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ABSTRACT 


Analysis of data from the Stanford University ELF and VLF experiment 
on 0G0-4 satellite has revealed new types of ion whistlers and also some 
new features of low to medium latitude ELF hiss not previously identified. 
The ion whistlers identified and interpreted are the l - proton whistler, 
the ion-cutoff whistler and combinations of both forms. The observations 
of ELF hiss permitted the identification of a hiss band and its interpre- 
tation in terms of propagation effects, using computer ray tracings. The 
band exhibits the previously observed low frequency cutoff characteristics 
of downward propagating hiss and also a peculiar upper cutoff in the 
vicinity of 600 Hz, which is nearly independent of the satellite altitude 
(~ 430 to 900 km) and latitude above ~ 15°. This hiss band is called 
band-limited ELF hiss, or simply BLH. It is observed between ~ 10° and 
~ 55° dipole latitude, where there is merging of the nearly constant 
upper cutoff and the increasing (with latitude) lower cutoff. Around 
~ 10° dipole latitude the BLH exhibits a sloping upper cutoff decreasing 
in frequency toward the equator (equatorial erosion). The BLH is seen 
most frequently from ~ 0600 to ~ 2200 LT, although spme examples have 
been found between 0200 and 0500 LT, with less intensity. In the range 

0 to 3 kHz and up to ~ 55° dipole latitude the BLH is the strongest 

-4 2/ 

signal observed. Its peak amplitude may reach 2 x 10 y /Hz between 

40° and 50° dipole latitude during daytime ( ~ 1000 LT). After prolonged 

periods of low magnetic activity (K p < 2) the upper cutoff may decrease 

to as low as 420 Hz but recovers about one day after a sharp increase in 

magnetic activity, such as a sudden commencement. The decrease and 

smearing of the lower cutoff of the background hiss, observed at latitudes 

greater than about 58°, indicates the presence of the light ion trough. 

- iii - SEL-71-055 



The similarity of the upper cutoff frequency on the two sides of the 

equator is taken as conclusive evidence that the source lies on the 

equator. An equatorial source near L = 4 and large initial wave normal 

angles were found adequate to explain the characteristics of the BLH as 

well as some of the observations of ELF hiss previously reported. Two 

possible generation mechanisms, Cerenkov and Doppler -shif ted cyclotron 

(DSC), are briefly examined. For the Cerenkov mechanism electrons of 

energy of the order of 1 kev are required, while for the Doppler-shifted 

cyclotron mechanism, much higher energies, of order 1 Mev are heeded. 

For a typical distribution of energetic particles the number density 

4 

of electrons for the Cerenkov mechanism is of order 10 times greater 
than that for the DSC mechanism. However, in either case a certain amount 
of emission coherency or amplification is necessary to obtain the observed 
power fluxes; 

The propagation of electromagnetic waves in a nonhomogeneous aniso- 
tropic medium is examined from the point of view of geometrical optics. 

In particular, the propagation of ELF waves in the magnetosphere is 
described in terms of the electron and ion densities and the intensity 
and inclination of the earth's magnetic field. The analysis of the 
variations of wave normal angle along the ray path is extended to include 
the effects of ions. A comparison of the relative importance of each of 
the above parameters in controlling the orientation of the wave normals 
is made in the region of the magnetosphere where most of the ion whistlers 
have been detected ; 
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I . INTRODUCTION 


A . PURPOSE 

The purpose of the present work is to discuss and interpret a 
certain class of electromagnetic waves occurring naturally in the earth's 
magnetosphere. They are found in the ELF range of frequencies (100 Hz 
to 3000 Hz) and are detected aboard artificial satellites. They include 
phenomena both whistlers and hiss and are interpreted in terms of the 
known properties of the propagation of electromagnetic waves in aniso- 
tropic dispersive media. 

The propagation of plane electromagnetic waves in a dispersive 
medium can be described in terms of the quantity known as refractive 
index. This quantity, defined by the ratio of the speed of light in a 
vacuum to the magnitude of the actual phase velocity of the wave con- 
sidered, is a function of the characteristics of the medium and, of 
course, of the frequency of the wave. When the medium is the earth's 
upper atmosphere, which is a composite semi-ionized gas, the propagation 
problem becomes particularly interesting due to the anisotropic charac- 
teristics of this medium. This means that the phase velocity of a plane 
wave is not only a function of position, but also of the direction of 
propagation . 

The anisotropy in the magnetosphere is caused by the presence of 
the earth's magnetic field, which puts some constraints on the movements 
of the charged particles (electrons and ions) that constitute the ionized 
atmosphere. When these particles are displaced from their positions of 
equilibrium, they will describe curved trajectories around magnetic 
field lines. The projection of this trajectory on a plane perpendicular 
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to the magnetic field direction is a circular motion with an angular 
frequency (gyrofrequency ) that depends on the particle's e/m ratio and 
on the field intensity. It is therefore to be expected that when these 
particles are set in motion by the electric field of a passing wave, 
they will have some effect upon the propagation of such a wave. When 
the frequencies of interest are much higher, than the ion gyrof requencies , 
the ions may be considered as stationary in comparison with the much 
lighter electrons, so their influence may be neglected. However, when 
the frequencies are in the range of the ion gyrof requencies , that 
approximation is no longer valid . The so called ELF range of frequencies 
falls under this category and the presence of the ions in the higher 
atmosphere (oxygen, helium and hydrogen) contributes appreciably to the 
dispersive properties of the medium. 

B. BACKGROUND OF THE PROBLEM 

The natural phenomena to be studied in this work belong to two 
broad categories known as: 

1. Whistlers, which are dispersed forms of electromagnetic energy 
produced in the atmosphere by lightning discharges. 

2. Hiss, which consists of continuous broadband noises believed to 
be produced by radiation from charged particles. 

In the case of whistlers, the energy released in the lower atmosphere 
by the sudden discharge must penetrate the atmosphere-ionosphere boundary. 
For the frequency range of interest in this work (100 to 3000 Hz), this 
boundary may be assumed to be sharp. For typical conditions of the lower 
ionosphere the decrease in phase velocity of the entering wave is of the 
order of 50:1. This decrease in phase velocity means, according to 
Snell’s law, that for any angle of incidence below the boundary, the 
waves will emerge on the other side with their wave normals practically 
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perpendicular to the boundary. The details of this penetration mechanism 
can be rigorously treated only by the application of full wave methods 
(see Budden [1961] and Scarabucci [1969]). However, the full wave 
approach has confirmed in general the simple geometrical optics reasoning 
used above and, unless one is particularly interested in studying the 
amplitude of the waves across the transition, the geometrical optics 
approach is adequate. For the present purposes, it is also satisfactory 
to consider the atmosphere -ionosphere boundary fixed at approximately 
100 km altitude. Until 1968, the only particular whistler traces whose 
characteristic dispersion could be attributed to specific ion effects 
were the proton and the helium whistlers (see Gurnett et al. [1965] and 
Barrington et al. [1966]). Since then, additional "ion whistlers" have 
been identified and interpreted by the author and their study is included 
in the present work. 

ELF and VLF hiss or continuous noises have been reported exten- 
sively in the magnetosphere, mainly from analysis of satellite data. 

The range of latitudes and frequencies is quite broad, from low latitudes 
to the polar region and from a few hundred Hz to dozens of kHz. In 
particular, the so-called auroral hiss, in the range of ~ 5 kHz to 
dozens of kHz and confined to the auroral zone has been treated exten- 
sively in the literature (see Jorgensen [1968]; Laaspere et al. [1971] 
and references in those papers) and is not the object of the present 
study. The kind of hiss of interest here (ELF hiss extending from a 
few hundred Hz to approximately 3 kHz) has been reported and discussed 

also extensively (see discussion in Chapter 4), but most investigators 

o 

have confined their attention to invariant latitudes greater than 35 
and/or high altitudes. In the low to medium latitudes and altitudes of 
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the order of 1000 km, the only particular feature of the ELF hiss 
established by previous investigators is the low frequency cutoff (see 
Gurnett and Burns [1968]). 

C. CONTRIBUTIONS OF THE PRESENT WORK 

The contributions of the present work can be stated briefly as 
follows: 

1. ' Expressions for the rate of variation of the wave normal 

direction along the ray path are developed and applied, 
including the effects of ions. This is an extension of previous 
work by Cerisier [1967] in which the effect of ions are not 
included. 

2. Certain new types of whistlers are identified and explained, 
including the transequatorial or l - proton whistler, the ion 
cutoff whistler and hybrid forms of both. 

3. New characteristics of an ELF hiss band (here called band 
limited ELF hiss or simply BLH) are found. Their interpre- 

. tation has led the author to suggest a model with the source 
located on the geomagnetic equator near L = 4. This new 
model provides alternative explanations of some previously 
reported properties of ELF hiss. 

D. PLAN OF THE WORK 

In Chapter 2, some introductory material is presented on the 
propagation of waves in anisotropic media, followed by the development 
and application of expressions for the rate of variation of the wave 
normal orientation along the ray paths . 

In Chapter 3, the 1 proton whistler, the ion cutoff whistler and 
hybrid forms of both phenomena are presented and interpreted using ray 
tracing techniques and concepts from Chapter 2. 

Chapter 4 deals with the observations and characteristics of ELF 
hiss and an interpretation also based on propagation properties. A 
range of wave normals at the source and a region of generation are then 
deduced and some comments are made about generation mechanisms based on 
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radiation from energetic electrons. 


Finally in Chapter 5 some conclusions are drawn and recommendations 
are made for future work. There are also four appendices where deriva- 
tions and auxiliary material are presented. 
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II. EFFECTS OF THE MEDIUM ON THE PROPAGATION OF ELF WAVES 


A . INTRODUCTION 


In the following analysis the medium is supposed to be a cold, 

collisionless, electrically neutral plasma composed of electrons plus 

three singly-charged positive ion species: oxygen, helium and hydrogen. 

Furthermore, this plasma is immersed in a static dipole type magnetic 

field. For the low frequency end of the spectrum, the refractive index 

p at a fixed point in such a med ium varies, as a function of frequency, 

as shown if Fig. 1 (for details see Stix [1962]; Smith and Brice [1964]; 

Gurnett et al. [1965]). The branches labeled R and L represent the 
2 

values of p corresponding to a wave normal direction along the static 

magnetic field, the R and L standing, respectively, for the right- 

hand and left-hand circularly polarized modes. The branches labeled e 

2 

(for extraordinary) represent the values of p for a wave normal 

direction perpendicular to the static magnetic field, with linear 

polarization (static and wave magnetic fields parallel). For intermediate 

2 

wave normal directions p lies inside the shaded regions. Using Stix's 
nomenclature and MRS units we can write: 



R . 
J 


X . 
J 


1 + s . Y . 
J 3 


L = 1 



L . 

J 


X . 
J 


1 - s . Y . 
3 3 


P = 1 - £ X . 

j J 


(0 


X . = 
J 


PJ 

2 

co 


where 


2 N j 6 

(0 . = — 

PJ m 


(plasma frequency) 


y = 

3 m . co 
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where 



where \jj is the angle between the static magnetic field and the wave 

& 

normal. From now on the angle will be simply called "wave normal 
angle . 

In Fig. 1 the frequencies labeled f to f are characteristic 

1 O 

frequencies and correspond to the various resonances (where p -* ®) 

and cutoffs (where p -* 0) . Frequencies f , f^ and f ? are, respectively, 

the gyrof requencies of the oxygen, helium and hydrogen ions. Frequencies 

f„ and f_ are called two-ion resonances and f is normally called 
2 5 o 

the lower hybrid resonance frequency (LHR). Similarly, frequencies f 

and fg are called two-ion cutoffs. Frequencies f^ and f x2 are 

the crossover frequencies, where the polarization of the R and L 

modes become linear (static and wave magnetic fields perpendicular) and 

interchange signs. The ion gyrof requencies depend only on the e/m ratios 

of each particle species and the intensity of the static magnetic field. 
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FIGURE 1. CURVES OF p 2 AS A FUNCTION OF FREQUENCY IN A 
MEDIUM COMPOSED OF ELECTRONS AND THREE ION 
SPECIES (OXYGEN, HELIUM AND HYDROGEN) . The 
several characteristic frequencies caused by the presence 
of these ions are indicated (f to f ). 

1 O 

For the models considered here their ratio is always 1:4:16. At a 
fixed point in space, frequencies f^, f 4 and f ? are constant and the 
others vary with the relative ion composition. It should be noted, 
however, that frequencies f and f are always situated between the 

fc O 

ion gyrof requencies f and f , while f and f are always situated 

1 4 O D 

between f 4 and Figures 2a-g illustrate what happens with the two 

ion resonances and cutoffs when the relative ion concentrations are 
varied . 

Due to the fact that the static magnetic field (and thus all the 
gyrof requencies) vary with the inverse cube of geocentric distance, an 
observer moving up or down inside this medium will see the group of 
frequencies f to f , respectively, decrease or increase monotonica lly . 

- 8 - 
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FREQUENCY 


3,4,5 6 


1 


8 


FIGURE 2a, b. CURVES OF y 2 AS A FUNCTION OF FREQUENCY FOR 
TWO DIFFERENT IONIC COMPOSITIONS. The several 
characteristic frequencies are indicated (f^ 
to fg)* In (a) the oxygen relative concentration is negligible 
when compared with the other two ions. In (b) the helium 
relative concentration is negligible and the hydrogen relative 
concentration is smaller than that of oxygen. 
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FIGURE 2c, d. IN (c) THE HELIUM RELATIVE CONCENTRATION IS 
NEGLIGIBLE AND THE HYDROGEN RELATIVE CON- 
CENTRATION IS LARGER THAN THAT OF OXYGEN. 

In (d) the hydrogen relative concentration is negligible. 
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4 


FREQUENCY 8 


1,2,3 


5,6,7 


FIGURE 2e,f. IN THESE EXAMPLES ONLY ONE OF THE ION SPECIES 

IS OF IMPORTANCE. In (e) the helium and oxygen 
relative concentrations are negligible. In 
(f) the hydrogen and oxygen are negligible. 
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FIGURE 2g. HERE THE ONLY ION OF IMPORTANCE IS THAT OF 
OXYGEN; THE HYDROGEN AND HELIUM RELATIVE 
CONCENTRATIONS ARE NEGLIGIBLE. 


2 

This means that the topology of the curves of p versus frequency for 

a fixed point in the magnetosphere is similar to the topology of curves 
2 

of p versus altitude for a constant frequency. Now, there will be 
characteristic levels in analogy with the former characteristic fre- 
quencies. Consequently, a point moving along the positive frequency 
axis in Figs. 1 or 2 may also be interpreted as a constant frequency 
wave moving upwards . 

At each point inside the magnetosphere one can define a vector 
with length equal to p and same direction as the wave normal. The 
tips of all possible refractive index vectors describe a surface that is 
called the refractive index surface. As the medium is uniaxial, the 
refractive index surfaces are figures of revolution around the static 
magnetic field direction. Figure 3 shows the intersection of a typical 
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FIGURE 3. THE RAY DIRECTION (S) IS ALWAYS PERPENDICULAR 
TO THE REFRACTIVE INDEX SURFACE. 

refractive index surface with a plane containing B. It can be shown 
(see Stix [1964]) that the direction of energy flow or ray direction S 
is normal to the refractive index surface at the point p where the 
vector p meets that surface. It can be also shown that: 

tana = -^|r (2 - 2) 

where a is the angle from the wave normal direction to the ray direction. 

>•4 ^ — * 

Therefore, at each point, the three vectors S, p and B are contained 
in the same plane. Furthermore, as dp/cty is positive for 0 < ijj < jt/ 2, 
the ray direction lies between the magnetic field and the wave normal 
directions . 

Finally, one may observe that each two-ion resonance divides the 
neighboring (p,f) space and its analogue, the (p,h) space, in two regions 
where the refractive index surfaces have different topologies. The region 
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immediately below f 2 or f 5 is characterized by closed refractive 

index surfaces, while the region immediately above f or f is 

2 6 

characterized by open refractive index surfaces. For large wave normal 
angles then, waves propagating in the two regions will have very different 
ray directions (see Fig. 4). Therefore, it can be expected that a wave 
starting its path at the base of the ionosphere between, for example, 
frequencies f ^ and f (see Fig. 1) and propagating upwards, will 
experience a noticeable change in ray path when it crosses the resonance 
level. This effect will be shown in Chapter 3 with the aid of ray tracings. 


B. FACTORS AFFECTING THE ORIENTATION OF THE WAVE NORMALS 

The propagation of ELF waves throughout the magnetosphere is 
affected by the presence of positive ions mainly through the cutoffs 
and resonances that they introduce in the dispersion relation. In 
order to understand the effects of the electron and ion densities and 
also of the earth's magnetic field on the orientation of the wave normals 
along the ray paths, the following expression, derived by Cerisier [1967] 
from the Haselgrove [1954] ray tracing equations will be used: 


djjt. 

ds 


cos a 




1 Su 
+ n Sb 


VB 


- vp 



(2.3) 


where 

ds represents an element of ray path 

Nj represents the electron and ion number densities 

B is the magnitude of the static magnetic field 

(3 is the inclination of the magnetic field relative to the 
dipole axis (see Fig. 5) 

k is a unit vector along the wave normal direction 

k is a unit vector perpendicular to the wave normal direction, 

X oriented along ijf + it/ 2 (see Fig. 5). 
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DIPOLE AXIS 



FIGURE 5. POLAR COORDINATE SYSTEM. All angles are 
measured clockwise from the references 
shown . 

The three terms inside the brackets are the contributions of the 
electron and ion densities, the magnitude of the dipole magnetic field 
and the inclination of this magnetic field. Equation (2.3) is derived 
in Appendix A. 

To apply Eq . (2.3) Cerisier used an approximate expression for the 
refractive index, neglecting the influence of the ions. In the following 
analysis the complete expression for p will be used. The derivation 
of the equations that follow can be found in Appendix A. 

The first term inside the brackets in Eq . (2,3), which corresponds 
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to the contribution of the electrons and ions can be written as: 


S - 


3n~ 

3 


VN 


s 

j 


1 d(l 

n 'Sr 


R . + 
J 


i Sn 

ii Sl j 




(2.4) 


where 




Sn. , Sn. 

J . . 1 3 t 

Sin 6 " r W COS 6 


(2.5) 


6 is the angle between the local vertical and the wave normal direction 
and (r,9) are the polar coordinates of a point in the magnetic meridian 
plane (see Fig. 5). 

Similarly, the second term in Eq . (2.3), which corresponds to the 
contribution of the magnitude of the static magnetic field (dipole), 
can be written: 


where 


Is VB 



CU b 2 ^ t 2 1 VB 
Sr R j ” c)l L j 1“ 



sin 6 - cos 


sin 0 cos S 
2 

1 + 3cos 9 


Finally, the third term, for the contribution of the inclination 
of the magnetic field is: 

2 

« 3 „ 1 + cos 9 

Vp • k = - sin 6 o 

1 + 3cos e 

The effectiveness of each of the three terms in Eq. (2.3) in determining 
the change in direction of the wave normals depends on the angle a 
between the ray and the wave normal. When a is zero, the effect is 
maximum and when a is it/2 the effect is zero. When a < «/2, cos a 
is positive and the resulting sign of di|//ds must come from the terms 
inside the brackets. If the wave normal angle \[i is measured from the 
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direction of B with the provision that jty| < «/ 2 , the sign will be 

positive (negative) when the wave normal k turns towards the inside 

(outside) of the field line. Because of the scalar products 7N • k 

j l 

appearing in the first term of Eq. (2.3), each of its four components 
will change signs whenever the wave normal crosses the direction of the 
respective gradient. The same is true of the second term when the wave 
normal crosses the direction of VB. As for the third term, expressed 
by Vp • k , the change in sign occurs around the vertical direction, 
because Vp is always horizontal for a dipole field. 

C. EVALUATION OF THE COMPONENTS OF diji/ds IN A MODEL MAGNETOSPHERE 
In this section, the three components of dijf/ds are evaluated in 
a model magnetosphere for a given ELF wave. The computation is made as 
a function of altitude, maintaining a constant wave normal direction in 
order to compare the magnitudes of the three terms in several regions of 
the (|_i,r) space around the characteristic levels of the hydrogen ion 
(where the wave frequency is equal to one of the characteristic frequencies). 
The density model used here is 1 composed of electrons and three ion species 
(hydrogen, helium and oxygen). The electrons and ions are assumed to be 
in diffusive equilibrium along magnetic field lines in the manner 
described by Angerami and Thomas [1964]. The temperature is the same 
for all particles and the parameters of the model are adjusted to obtain 
a reasonably good fit with an electron density profile from Alouette II 
in the fall of 1967. It should be stressed at this point that no claims 
whatsoever are made here about how representative of actual conditions 
those parameters are. The D.E. model is chosen mainly because it is well 
behaved mathematically. Details of the model and derivation of expressions 
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for ^ and appearing in Eq. (2i5) can be found in Appendix B. 

Figure 6 shows the electron density profile of the model between 
450 and 2000 km of altitude. The dots represent values of density 
measured by Alouette II in August 5, 1967 at approximately 2000 LT . The 
parameters of the D.E. equations are also shown in the figure and reference 
to Appendix B should be made for the meaning. The range of altitudes of 
450 to 2000 km covers the regions where most of the so-called ion whistlers 
have been detected . 



FIGURE 6. ELECTRON AND H + DENSITY PROFILES FOR THE DAYTIME 
IONOSPHERE MODEL USED THROUGHOUT THE PRESENT WORK. 
The parameters indicated are defined in Appendix B. 



The dispersion relation admits two solutions, or modes, that may 
or may not correspond to propagating waves in the medium. Whenever the 
two modes can propagate simultaneously, the terminology "slow" and "fast.' 
can be applied unambiguously in correspondence with the respective values 
of phase velocity or refractive index. However, when one of the modes 
is evanescent, the same terminology will be extended here as shown in 
Fig. 7 around the hydrogen characteristic frequencies. In this context:, 
a typical representative of the slow mode is the proton whistler and a 
typical representative of the fast mode is the ion-cutoff whistler to be 
studied in the next chapter. 

( 

2 I 

U 



FIGURE 7. TERMINOLOGY USED FOR THE TWO PROPAGATING MODES 
IN THIS WORK. 
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Figures 8 to 11 show curves of the three components of d^/ds . 
between 450 and 2000 km and a latitude of 20°, for a wave frequency of 
300 Hz and a constant wave normal angle of -10°. Because a logarithmic 
scale is used, the magnitude only of each term is plotted and the signs 
are indicated along the curves. Figures 8 and 9 correspond to the slow 
mode and extend up to 1150 km altitude, very close to the resonance level 
where |_i =* «. Above this altitude the mode is evanescent. Above 450 km 
the particle density term is positive and is dominated by the electron 
density contribution. In the vicinity of the crossover altitude, the term 



RAD/km 


FIGURE 8. CURVES OF THE COMPONENT OF dijj/ds DUE TO 
THE ELECTRON AND ION DENSITIES (N ) , FOR A 
WAVE PROPAGATING IN THE SLOW MODE . 
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FIGURE 9. CURVES OF THE COMPONENTS OF dty/ds DUE TO 
THE MAGNITUDE (B) AND INCLINATION (p ) OF 
THE EARTH'S MAGNETIC FIELD, FOR A WAVE 
PROPAGATING IN THE SLOW MODE. 


decreases very fast and changes sign. From this point upwards, the 

dominant contribution is that of the hydrogen ions. If the wave normal 

angle were zero, the resonance would occur exactly at the proton gyro- 

frequency altitude, which is independent of density variations. Therefore, 

if i|f = 0 the particle density term would remain finite up to and through 

the resonance altitude. However, if the wave normal angle is different 

1 

from zero, it can be shown that — ■=» °o at resonance, which makes the 

particle density term also infinite. For small angles (such as in the 


examples) the magnitude of 
SEL-71-055 . 
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decreases with increasing altitude up to the immediate vicinity of the 
resonance artd then grows very fast to infinity, as indicated in Fig. 8. 

The terms corresponding to the magnitude and inclination of the magnetic 
fields are illustrated in Fig. 9. The contribution of the magnitude of 

■*4 

B is negligible in comparison with the particle density term up to the 
neighborhood of the crossover altitude and then becomes entirely dominant. 
However, unlike the particle density term, this one does not change sign 
anywhere and its magnitude always grows to infinity at resonance, inde- 
pendently of \|f . The term due to the inclination of B is independent 
of p arid decreases steadily as i/r; 

Corresponding curves for the fast mode are shown in Figs. 10 and li. 
The fast mode exists only above the two-ion cutoff altitude where the 
refractive index goes to zero* making the particle density term and also 
the term due to the magnitude of B tend to infinity. Between the cut- 
off altitude and approximately 1005 km both terms have the same order of 
magnitude and same sign. In this region, the particle density term is 
dominated by the contribution of the hydrogen ions. The term due to 
the magnitude of B continues to decrease untii about 1020 km but 
without changing sign; Above that altitude the variations are much 
less pronounced and its magnitude is comparable to that of the third 
term (V|3 • k) . The particle density term is the most important between 
~ 1020 and 2000 km but decreases faster than the other two with increasing 
altitude. Therefore, one can expect this term to become negligible at 
larger altitudes. 
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FIGURE 10. CURVES OF THE COMPONENT OF diji/ds 

DUE TO THE ELECTRON AND ION DENSITIES 
(N.), FOR A WAVE PROPAGATING IN THE 
FAliT MODE. 
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III. ICN WHISTLERS 


A . INTRODUCTION 

The natural whistler phenomena to be examined in this chapter may 
all be loosely called "ion -whistlers , " meaning whistlers whose spectral 
shape is mainly a result of ion effects on the dispersion relation. 

These whistlers have been observed in data from the Stanford University 
VLF experiment flown on the 0G0-4 satellite. This is a quasi-polar 
orbiting satellite ( ~ 86° inclination) with perigee of ~ 400 km and 
apogee of ~ 900 km (orbital period ~ 98 minutes). The signals are 
detected in a magnetic loop antenna which feeds a number of receivers 
serving different purposes. Of interest here are the two broadband 
receivers used for spectrographic analysis. The so-called ELF receiver 
covers the band of ~ 15 to 300 Hz with a gradual rolloff between 300 
and 1000 Hz. The VLF receiver covers the band of 300 Hz to 12.5 kHz. 

The signals are telemetered to ground and recorded in magnetic tape at 
the NASA telemetry stations and later they are processed at Stanford 
using a Rayspan real time spectrum analyzer (see Helliwell et al . [1961]). 
The final results are presented as frequency versus time spectra in 
strips of 35mm film with the time running alongside the film. 

Following the conventional terminology used for whistlers detected 
in satellites (see Smiths and Angerami [1968]), the events to be reported 
here will be called 0 + (zero plus), 1 (one minus), 1 + (one plus), etc. 
according to the relative positions of the satellite and the source with 
regard to the magnetic equator and also according to the direction of 
propagation at the time of detection. 
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B. THE 1 PROTON WHISTLER 

This peculiar trace has been observed in some equatorial passes 
of 0G0 4 [Muzzio, 1969] and a few examples are indicated by arrows in 
the spectrograms of Fig. 12. The two strips represent sections of the 
same pass detected at the Santiago telemetry station. The frequency 
range is 0 to 1.25 kHz and the time scale is indicated at the upper 
part of each strip. The traces are very similar to, and appear among 
regular or 0 + proton whistlers from which they differ mainly by the lower 

value of their frequency of infinite delay time (f ) . In addition, when 

00 

the satellite moves away from the magnetic equator, the value of f 

CO 

for the events labeled 1 decreases. The following analysis will show 
that the traces are proton whistlers propagating from the opposite 
hemisphere approximately along field lines. The name 1 proton whistler 
is therefore appropriate. The decrease in f can be appreciated in 

CO 

the top panel of Fig. 13 until ~ 0736, when the 1 type traces disappear. 
Near 0735 some 0 + proton whistlers from the southern hemisphere begin to 
be visible. Sections A and B of the top spectrogram are shown expanded 
in time at the middle and lower left, where individual events can be 
resolved. The approximate ray paths of events 1 (0 + ) and 2 (1 ) are 
indicated in the sketch at right. In the 1 type traces, frequencies 
above ~ 150 Hz start their path at the base of the ionosphere already 
above the helium gyrofrequency level (frequency in Fig. 2e in 

Chapter 2). These waves are propagating in a region of closed refractive 
index surfaces. The ray direction is situated between the wave normal, 
which is initially vertical, and the direction of the static magnetic 
field. As the waves progress upward, they will approach the two-ion 
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FREQUENCY (kHz) 


SNT 0G0-4 19 SEP. 1967 (southbound) 

0652 (-15° Dip. Lot.) 


1.25 


0 


1.25 


0 


FIGURE 



12. EXAMPLES OF 1 PROTON WHISTLERS. Notice the 

difference in values of f between the 0 + 
_ 00 
and 1 events. 
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SNT 0G0-4 17 SEP 1967 

0731 0732 0733 0734 0735 0736 0737 UT 
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FIGURE 13. ILLUSTRATION OF THE DECREASE IN VALUE OF OF THE 1 PROTON WHISTLERS 

(LABELED 2 IN THIS FIGURE) WITH INCREASING LATITUDE. The approximate 
ray paths of the 0 + (labeled 1) and 1~ events are indicated in the 
sketches . 


resonance level (f,_), above which the refractive index surfaces open 
and the guiding effect of the static magnetic field becomes more and 
more important. To illustrate this change in ray path direction across 
the two-ion resonance level the ray tracing results of Fig. 14 are very 
useful. The frequency used is 200 Hz and the ray tracing starts at 
100 km altitude and 19° dipole latitude. The ionosphere model is the 
same as used in Chapter 2. The wave normal is initially vertical and 
the arrows along the ray path indicate subsequent orientations. The 
transition from closed to open refractive index surface is also indicated 
and the corresponding change in ray path direction can be noticed. When 
the ray crosses the equatorial plane the alignment with the field is 
already very good. The variations of wave normal angle along the ray 
path can be roughly compared with the values of dijj/ds for the slow 
mode plotted in Figs. 8 and 9. From the values of wave normal angle 
marked along the ray tracing of Fig. 14, one can observe the change in 
sign of dty/ds close to the crossover level. 

The upper frequency limit for a wave propagating in the slow mode, 
or similarly, the upper altitude limit for a fixed frequency, is given 
by the resonance condition: 

tan^tli = - ^ (3.1) 

where : 

S = i (R + L) (3.2) 

When the wave normal angle varies from 0° to 90° (see Stix [1962]) 
the resonance level decreases from the' proton gyrof requency to the two- 
ion resonance level. However, if j ijj | is less than about 75 the 
resonance remains close to the proton gyrof requency level. Consequently, 
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for waves in the slow mode with ray paths similar to the one in Fig. 14, 
the upper frequency limit across the equatorial plane will be close to 
the minimum proton gryof requency along the ray path, instead of the local 
value at the satellite. This will be the equatorial value of gyro- 
frequency along the field line passing through the satellite. Therefore, 
as the satellite moves away from the equator the minimum gyrof requency 
along the ray path decreases, which is precisely the picture in Fig. 13. 
For each satellite altitude, there is a value of latitude beyond which it 
is impossible to detect 1 proton whistlers. This limit can be under- 
stood from Fig. 15, which represents curves of proton gyrof requency , 
crossover frequency and two-ion resonance frequency along a field aligned 
ray path containing the equatorial region. The characteristic frequencies 
of the other two ions have been omitted, for simplicity, but they could 
be used to explain the eventual occurrence of similar whistler traces 
at lower frequencies. The vertical lines marked (I) and (2) represent 
two satellite positions along the same field line and the corresponding 
spectral forms received at these two positions are illustrated in the 
two spectrogram sketches in the same figure. The minimum values of the 
characteristic frequencies occur at the equator, which is the highest 
point along the ray path. At higher latitudes the ray path is lower in 
the ionosphere. The corresponding decrease in hydrogen concentra tion 
makes the crossover and two-ion resonance frequencies approach the proton 
gyrof requency , enhancing the probability of coupling between modes. 

This region of strong coupling is indicated in the figure. In the 
following analysis, waves in the slow mode are supposed to be propagating 
from left to right in the diagram and the effect of heavier ions will be 
disregarded. For simplicity, coupling near the crossover frequency will 
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FREQUENCY 





FIGURE 15. CURVES OF THE CHARACTERISTIC FREQUENCIES ALONG 
A FIELD-ALIGNED TRANS EQUATORIAL RAY PATH AND 
CORRESPOND ING SPECTRA OF WO 1" PROTON WHISTLERS 
RECEIVED AT POSITIONS (l) AND ( 2 ) . 
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be disregarded also, except in the shaded regions where the characteristic 
frequencies are very near each other. With these restrictions, only 
frequencies below f will be able to reach either one of satellite 
positions (T) or (g) in the slow mode. Depending on the satellite 
position however, there will be another limit on the low frequency side. 

To understand this limit, Figs. 14 and 15 have to be examined simultaneously. 
After crossing the equatorial plane the rays move downward and the wave 
normal angle approaches -90°. In these circumstances, when the two-ion 
resonance level is attained the waves will be reflected back along a 
very near ray path. This mechanism is illustrated in Fig. 16 by the use 
of a "Poeverlein construction." On the other hand, the energy of the 
waves will probably be damped by lossless temperature effects across 
this resonance (see Stix [1962], Chapt. 3). For 1 proton whistlers 
therefore, the minimum frequency observable will be either the minimum 
crossover frequency along the ray path (satellite position (l) in Fig. 15), 
or a frequency very near the local two-ion resonance frequency (satellite 
position (5) ) , whichever is bigger. This constitutes a "window" limiting 
the frequency range of observable 1 proton whistlers . As the satellite 
moves to higher latitudes this window becomes narrower to a point where 
it collapses. This limiting latitude is indicated in Fig. 15 and it 
happens where the local two-ion resonance frequency equals the equatorial 
value of proton gyrof requency . Above that latitude, a satellite in 
circular or near circular orbit will not be able to detect any additional 
1 proton whistlers. It should be understood that the results indicated 
by Fig. 15 are only approximate, since the damping due to temperature 
effects near the resonances have not been taken into account quantitatively. 
The above mechanism is very adequate to explain all the features of 1 
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FIGURE 16. MECHANISM OF REFLECTION OF A WAVE IN THE SLOW 
MODE WITH LARGE WAVE NORMAL ANGLES, ACROSS THE 
TWO- ION RESONANCE LEVEL. 


proton whistlers illustrated in Fig. 13, including their disappearance 
o 

around -19 dipole latitude. 

Finally, one comment about the polarization of these whistlers. 

Due to the fact that the crossover level marks a region of polarization 
reversal, the actual polarization of the received 1 proton whistlers 
will depend on the relative position of the satellite among the 
characteristic frequency curves. For example, in Fig. 15, the 1 proton 
whistler received at satellite position (T) is almost entirely left-hand 
polarized (except for the narrow frequency range below the local cross- 
over frequency), but the one received at position (5) is entirely right-hand 
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C. THE ION CUTOFF WHISTLER 

The previous section dealt with effects of ions upon waves in the 
"slow mode." The ion cutoff whistler, on the other hand, is the result 
of ion effects on the "fast mode." It was discovered in recordings of 
the 0G0-2 and OGO-4 satellites and interpreted on the basis of a reflection 
of downward propagating waves close to the two-ion cutoff frequency 
[Muzzio, 1968a,b], hence the denomination ion cutoff whistler. Some 
examples of this event are shown in the spectrograms of Fig. 17. Events 
(a) and (b) were detected over Australia, around a dipole latitude of 
-40° and events (c) and (d) were detected at Santiago (Chile), at a 
dipole latitude of -21°. The measured frequency -time curve for these 
whistlers fits a constant dispersion (Eckersley type) curve over most of 
the range 0-10 kHz. However, at the low frequency end the trace deviates 
towards increasing time delays, reaches a turning point or minimum 
frequency and then rises in frequency to some extent. Sometimes the 
trace is seen down to the minimum frequency only, the return portion 
being absent. 

Comparison of the observed whistlers with Eckersley curves indicate 

l/2 

dispersions of about 40 sec for the examples of Fig. 17a, b, and about 

l/2 

30 sec for the examples in Fig. 17c, d. These figures suggest the 

origin to be in the conjugate hemisphere, in other words, they are 1 

whistlers. The values of proton gyrof requency at the two locations shown 

in Fig. 17 are: 570 Hz for (a) and (b), and 300 Hz for (c) and (d). 

The corresponding minimum frequencies of the whistlers are lower than 

the above values of gyrof requency , as can be seen in Fig. 17. The rapid 

decrease in the refractive index of the fast mode when approaching the 
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FIGURE 17. EXAMPLES OF ION CUTOFF WHISTLERS DETECTED 
AT TWO DIFFERENT LATITUDES. In the upper 
and lower panels, f and F are, 
respectively, the minimum frequency and tlie local proton 
gyrof requency . 

two-ion cutoff level creates a reflection condition for waves propagating 
downwards. This reflection can be understood by the application of 
Snell's law or by the Poeverlein construction illustrated in Fig. 18. 

As a matter of fact, the analysis of the rate of change of wave normal 
angle along the ray path undertaken in Chapter 2 shows a very steep 
increase in d\j;/ds when the two-ion cutoff level is approached from 
above. The peculiar low frequency end of the whistlers in question 
suggests therefore, that the rays are being reflected somewhere below 
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POEVERLEIN CONSTRUCTION 


FIGURE 18 . REFLECTION OF A WAVE IN THE FAST MODE 

(PROPAGATING DOWNWARD) NEAR THE TWO-ION 
CUTOFF LEVEL. 


the satellite and cross the satellite orbit once more on the way up. 

Figure 19 illustrates the above argument. The curves marked f and 

c 

F represent, respectively, the two-ion cutoff and the proton gyro- 
H 

frequencies as a function of altitude. In the absence of coupling 
around the crossover level, a wave of frequency f propagating downward 
in the fast mode passes by the satellite at h^, is reflected near the 
cutoff level and may be detected a second time on the way up. The time 
interval At between the reception of the downward and upward signals 
depends on the distance Ah between the satellite altitude and the 
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reflection altitude, and also on the group velocity. The minimum 

frequency observed (f ) must be reflected at the satellite altitude. 

m 

Frequencies lower than f are reflected above the satellite and can 

m 

not be detected. This mechanism of reflection of ion cutoff whistlers 

proposed by the author in 1968 has been experimentally verified later 

by other investigators (see Mosier [1969], Fig. 5). The maximum frequency 

observed in the "tail" (f ), must correspond to the region where f 

M c 

and F come so close together that coupling occurs for any value of 
H 

wave normal angle (on the subject of coupling, see Jones [1969a, bj). 

Figure 19 also illustrates what happens when the satellite altitude is 
decreased down to the region of strong coupling. The tail of the ion 
cutoff whistler is progressively shortened until it can not be observed 
anymore. An example of this effect can be seen in the spectrogram of 
Fig. 20, where the tails of the ion cutoff whistlers are barely visible. 
Note the low altitude ( ~ 430 km) in contrast with both examples of 
Fig. 17 where the satellite altitude was above 850 km. 

One other effect that may limit the observability of the low 
frequency end of ion cutoff whistlers is the attenuation due to Coulomb 
collisions in the ionized medium. When collisions are involved, the 
refractive index becomes a complex quantity, a fact that is well known 
and will not be discussed in detail here (see for example Budden [1961]; 
Ratcliffe [1962]). It is possible to write a differential equation for 
the attenuation due to the imaginary part of the index and to incorporate 
this equation into a ray tracing program, for example. Following Budden 
[1961] and Kimura [1966] the rate of attenuation along the whistler ray 
path can be written as: 

~ = 181.92 f ^ 

dt p 
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FIGURE 20. EXAMPLES OF ION CUTOFF WHISTLERS WHERE 

THE RETURN PART OF THE TRACE IS NOT SEEN. 
Notice the low satellite altitude. 


where 

A is the attenuation in dB 

f is the frequency in Hz . 

X is the imaginary part of the refractive index 

p is the real part of the refractive index 

t is the time of phase travel along the ray path in sec. 

In the absence of collisions, the refractive index (real) drops 
to zero at the cutoff frequency. When collisions are present and the 
index becomes complex, neither the real nor the imaginary parts can be 
zero anymore and in reality there are no more cutoffs. However, if the 
collision frequencies are small enough, p may decrease to a sufficiently 
low value, close to the former cutoff, so that reflection of waves can 
take place. At the same time x increases, so the resulting attenuation 
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may attain very large values even before reflection occurs. We are 
interested here in momentum transfer collisions between charged particles 
(electrons and ions). Expressions for the collision frequencies can be 
obtained from Nicolet [1953]: 
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V, 


represent the collision frequency of the electrons 
with the three ion species considered throughout this 
work . 


v i 


represent the collision frequency of each ion species 
with the other two. Collision of ions with electrons 
are disregarded because only a small fraction of the 
ion momentum is transf ered . 


The above expressions were tested in a ray tracing program with a simulated 
ion cutoff whistler and the resulting attenuations versus reflection 
altitude are plotted in Fig. 21. The attenuations are measured from the 
equatorial plane and only values above 1 dB are plotted. The ionospheric 
model is the same daytime model used in Chapter 2. For frequencies that 
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FIGURE 21. ATTENUATION VS ALTITUDE FOR WAVES IN THE FAST MODE PROPAGATING FROM THE EQUATORIAL 
PLANE. The attenuation is produced by the effect of Coulomb collisions. The 
initial wave normals are horizontal and the magnetic field intensities are 
applicable to a region near Santiago. 



are reflected in the vicinity of 800 km altitude (or less), the attenu- 
ation is considerable, even close to the crossover level. Therefore, 
if the return trace is not seen in an ion cutoff whistler, one can not 
be sure if the apparent cutoff is caused by a reflection mechanism or 
by attenuation due to losses in the medium. There is still another 
possibility for the apparent cutoff which has to do with coupling condi- 
tions in the upward part of the ray path in the conjugate hemisphere. 

This condition will be dealt with in the next section regarding hybrid 
forms . 

D. HYBRID FORMS 

Having examined the propagation characteristics of the slow and 
fast modes separately in the last two sections, one is naturally led to 
ponder about the possibilities of simultaneously observing both modes 
associated with the same 1 whistler trace. In fact, these hybrid forms 
do appear in some satellite records, and an example is shown in Fig. 22. 
It represents a sequence of spectra from a single pass of OGO 4. The 
first two panels show only ion-cutoff whistlers, but as the vehicle 
approaches the magnetic equator, some 1 proton whistlers begin to 
appear, vestigially in (c) but stronger in (d). In (a) and (b) the 
satellite is too far from the equator to receive any 1 proton whistlers 
(see explanation in Section B) . 

The spectral shape of these hybrid forms can be understood very 
easily by using plots of the characteristic frequencies along ray paths, 
as was done previously. In Fig. 23 the curves marked F^, f^ and f 
are, respectively, the proton gyrof requency , crossover and two-ion cutoff 
frequencies along a transequator ial ray path. In reality, the ray paths 
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FREQUENCY (KHz) 


SNT 0G0-4 5 AUG. 1967 

22I8 = 20UT 2219=20 UT 



2221 = 30 UT 2222 = 40 UT 



ALTITUDE > 870 km 


FIGURE 22. TRANSITION FROM SIMPLE ION CUTOFF 

WHISTLERS (IN (a) AND (b)) TO HYBRID 
FORMS (IN (c) AND (d)) WHEN THE 
VEHICLE APPROACHES THE EQUATOR. 
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FIGURE 23. CHARACTERISTIC FREQUENCIES ALONG A TRANSEQUATORIAL RAY PATH AND 
SPECTROGRAM OF THE CORRESPONDING l” HYBRID FORM. 



for the two modes are not the same because of the very different shapes 
of the respective refractive index surfaces. However, this fact is not 
very important for the present discussion and would require eight curves 
in the diagram, an unnecessary complication. Therefore, the curves in 
Fig. 23 represent the characteristic frequencies along a hypothetical 

ray path, which is the same, for both modes. The regions bounded by 

/ 

these curves and labeled by the encircled numbers exhibit the following 
properties : 

1. Region (l) is forbidden to the fast mode and region (5) is 
forbidden to the slow mode. 

2. In the upward path, region ( 2 ) is inaccessible to the fast 

mode and region (5) is only accessible if there is coupling 

between modes across f . 

x 

3. In the downward path, the transition from (3) to ( 2 ) is allowed 
to the fast mode and there will be reflection in the vicinity 
of f c . There is also the possibility of coupling across f^. 

4. In region ( 2 ) the fast mode is left-hand polarized and the 
slow mode is right-hand polarized. Conversely, in region ( 3 ) 
the fast mode is right-hand polarized and the slow mode is 
left-hand polarized. Therefore, when there is coupling, or 
changing of mode, the polarization is conserved. 

5. The fast mode in region ( 2 ) and the slow mode in region © 
have comparatively low values of group velocity due to the 
large values of dp/df. 

Using the above properties in conjunction with Fig. 23, the following 
interpretation can be given to the hybrid forms: Frequencies below the 

minimum proton gyrofrequency along the path, while in the slow mode, 
will behave exactly as discussed in Section B (1 proton whistlers). 

Above a certain frequency indicated in the figure, mode coupling becomes 
possible around f in the upward path and so the fast mode can be 
excited. From this point on, the fast mode will propagate as an ion 
cutoff whistler and everything that was said in the previous section is 
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applicable here. The probability of mode coupling around f increases 

x 

when the characteristic curves approach each other, e.g., for higher 
frequencies (see Jones [1969a]). Ignoring for the moment the possibility 
of coupling in the downward path, and assuming that there is coupling in 
the upward path for frequencies equal or larger than the two-ion cutoff 
frequency at the satellite, the corresponding whistler trace received 
may exhibit the typical hook form (subjected to the problem of attenuation 
examined in the previous section). However, if conditions are such that 
coupling only occurs at higher frequencies, then the fast mode trace may 
not show any appreciable deviation from an Eckersley type dispersion 
curve and the cutoff seen will have nothing to do with a reflection 
phenomenon. These two situations are illustrated in the sketches of 

c 

Fig. 24 where, in the left are seen typical "coupled” and "uncoupled" 
electron-proton whistler pairs detected in the upward or 0 + path, and 
on the right, the corresponding forms after crossing the magnetic equator. 
In the uncoupled example the cutoff seen in the fast mode across the 
equator is the same as in the 0 + trace and is not caused by reflection. 

The frequencies in the "gap" stayed in the slow mode and became trapped 
in region (D (see Fig . 23) . 

If coupling in the downward path is taken into account, the energy 
could go from the fast to the slow mode across f before reflection 
close to the two-ion cutoff and a spectral form approximately equal to 
the sketch in Fig. 25 would result. To the author's knowledge, however, 
this kind of trace has not been observed yet. 
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FREQUENCY FREQUENCY 


COUPLED FORM 



TIME TIME 


UNCOUPLED FORM 



TIME TIME 


FIGURE 24. COUPLED AND UNCOUPLED ELECTRON -PROTON WHISTLER 
PAIRS SEEN AS 0 + AND 1" EVENTS. IN THE LOWER 
PANELS, f . INDICATES THE FREQUENCY ABOVE 
WHICH COUPtlNG STARTS FROM THE SLOW MODE TO 
THE FAST MODE IN THE UPWARD PATH. 
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FIGURE 25. EXAMPLE OF COUPLING FROM THE FAST MODE 
TO THE SLOW MODE IN THE DOWNWARD PATH. 


SEL-71-055 


50 


IV. ELF HISS 


A. INTRODUCTION 

ELF hiss or continuous noise in the range of a few hundred Hz up 
to about 2 or 3 kHz has been observed in the magnetosphere and reported 
by a number of authors (Guthart et al. [1968]; Taylor and Gurnett [1968]; 
Gurnett and Burns [1968]; Russell et al. [1969]; Dunckel and Helliwell 
[1969]; Mosier [1971]). Their work, however, deals almost exclusively 
with high latitude and/or high latitude events. 

This chapter presents and discusses observations of ELF hiss from 
the Stanford University VLF experiment flown on the 0G0-4 satellite 
( ~ 430 up to 900 km altitude) and covering the entire range between 
±84° invariant latitudes. Analysis of real time spectrogram data 
acquired at the Rosman (North Carolina), Quito and Santiago (South 
America) telemetry stations extending from September 1967 to March 1968 
revealed the occurrence of a hiss band with characteristics not previously 
identified. The band is seen from low to medium latitudes and will be 
referred to as band limited ELF hiss or simply BLH. Examination of 
several records from the telemetry stations at Johannesburg (South 
Africa) and Orroral (Australia) also revealed the same kind of hiss 
but a more thorough analysis of data from these stations has yet to be 
conducted . 

Furthermore, amplitude information from the digital stepping sweeping 
receivers of the same experiment was gathered from 14 satellite passes 
extending the range to ±84° invariant latitude. However, the results and 
conclusions from this sample are to be considered as preliminary until a 
more extensive study can be undertaken. 
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Based on wave propagation properties, it is proposed that the BLH 
is generated at large wave normal angles in the equatorial region near 
L = 4. This model can be used to explain the characteristics of the. 

BLH and also to interpret previous observations of ELF hiss reported by 
Russell et al. [1969] and Mosier [1971]. Two mechanisms for the generation 
of the BLH based on radiation from energetic electrons are considered. 
However, in order to explain the observed intensities, some form of wave 
growth or a coherency mechanism is required. 


B. OBSERVATIONS 

The ELF hiss detected by 0G0 4 exhibits the typical low frequency 
cutoff described by Gurnett and Burns [1968] and attributed to the two- 
ion cutoff in the dispersion relation. It also shows the following 
remarkable characteristics illustrated in the spectrograms of Fig. 26: 

1. The presence of a high intensity hiss band at low to medium 
latitudes below ~ 600 Hz (BLH) . When spectrograms covering 
a broad latitude range are examined, one generally observes 
ELF hiss extending from the lower cutoff (associated with 
the two-ion cutoff frequency) up to 2 or 3 kHz and from the 
equatorial through the auroral region (auroral hiss). In such 
a spectrogram, however, the BLH stands out as the strongest 
and steadiest signal in its latitude range. The hiss above 
the BLH is irregular and with a variable spectrum from day to 
day. Figures 26 and 28 illustrate different aspects of the 
background hiss and the BLH. 

2. The upper cutoff of the BLH in the vicinity of 600 Hz, which 
seems to be independent of latitude and is virtually the same 
in both hemispheres . 

3. The "equatorial erosion" or sloping upper cutoff decreasing in 
frequency toward the equator, near ±10° dipole latitude. It 
affects both the BLH and the background hiss. An equatorial 
erosion with a similar spectral cutoff has been previously 
observed in whistlers [Scarabucci, 1970] but, according to 
that author, it is produced mainly by absorption in the lower 
ionosphere. There are examples of simultaneous occurrences 

of erosion in the BLH and in whistlers, although the frequency 
ranges are different. Also, in some instances, the BLH exhibits 
the erosion and the whistlers do not. 
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Prior to this work, observations of ELF hiss by other investigators 
had established only that the lower frequency was controlled by the two- 
ion cutoff [Gurnett and Burns, 1968]. Reported observations were con- 
fined to invariant latitudes greater than ~ 35°, which precluded most 
investigators from observing the equatorial erosion. To the author’s 
knowledge, there are no references to the virtually constant upper 
cutoff of the BLH with latitude. 


The upper and lower cutoffs of the BLH taken from the spectrograms 
shown in Fig. 26 are plotted vs the absolute value of dipole latitude 


in Fig. 27. Although the lower cutoff is markedly lower in the southern 
hemisphere (due to the lower values of local magnetic field near the 
South Atlantic anomaly) the upper cutoff remains essentially the same in 


both hemispheres 
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FIGURE 27, 


VALUES OF THE UPPER AND LOWER CUTOFFS OF THE BLH 
SHOWN IN THE SPECTROGRAMS OF FIG. 26, PLOTTED VS 
DIPOLE LATITUDE. Although the lower cutoffs at 
conjugate latitutdes are markedly different due mainly to 
differences in local magnetic field, the upper cutoffs are 
essentially the same. 
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Figure 28 depicts two other spectrogram similar to the ones in 


Fig. 26 but recorded at higher latitudes by the Rosman telemetry station. 

Hiss and chorus up to ~ 2 kHz beyond 65° invariant latitude can be seen. 

On this record, the BLH is observed up to ~ 55° invariant latitude, 

where the upper and lower cutoffs merge. Above this latitude the lower 

cutoff can still be seen in the background hiss until its gradual dis- 

o 

appearance around 58 invariant latitude. This disappearance of the 
lower cutoff is consistent with the reported decrease in the H + relative 
concentration at high latitudes (see Taylor et al. [1969]), which greatly 
increases the probability of coupling for downgoing waves around the 
crossover altitude (for a study of the variations of the coupling co- 
efficient with ion concentration, see Gurnett et al . [1965); Jones [1969a, b) ; 
Rodriguez and Gurnett [1971]). To the author's knowledge this phenomenon 
has not been identified previously, although it can be seen in Fig. 5 of 
Rodriguez and Gurnett [1971]. The light ion trough seems to be a mani- 
festation of the plasmapause . In this case, the decrease and smearing 
of the lower cutoff. of the ELF hiss can provide a convenient way of 
determining the plasmapause location. 

In what follows we will be discussing mostly the properties and 
characteristics of the BLH. Examination of 270 occurrences (from spectro- 
grams of the OGO-4 broadband receivers) of hiss between September 1967 
and March 1968 from Rosman, Quito and Santiago revealed the following 
characteristics ; 

1. The BLH was seen in nearly all passes between 0600 and 2200 LT. 

No occurrences of the BLH were observed in the broadband data 
after ~ 2300 LT, although a few examples were found in the 
digital data (with smaller amplitude). The smaller amplitude 
of the nighttime events may be the reason why they have not 
been observed in the broadband receivers. A more extensive 
examination of the digital amplitude data is, therefore, recommended. 
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JAN 1968 



The gradual disappearance of the lower cutoff of the hiss between 58 a 
is probably due to the H + trough reportedly observed at the plasmapause and could 
provide a convenient way to determine its location at low altitudes. 


2 . 


The equatorial erosion is present in all records that cover 
the region around ±10° dipole latitude, and there is a marked 
symmetry in the phenomenon with respect to the magnetic equator. 

3. The upper cutoff of the BLH in the vicinity of 600 Hz is nearly 
independent of satellite latitude above ~ 15° and its sharp- 
'ness shows variations from pass to pass. However in ~ 40% of 
the passes examined the sharpness is comparable to that in 
Figs. 26 and 28. After prolonged periods of low magnetic 
activity the upper cutoff frequency decreases to values as 
low as 420 Hz, but recovers about one day after a sharp increase . 
in activity (such as a sudden commencement) . Figure 29 illus- 
trates the change in upper cutoff after the quiet period from 
7 to 9 March 1968. The spectrogram on the left (March 2) is 
typical of the period indicated from 1 to 6 March and has an 
upper cutoff of ~ 590 Hz. The spectrogram in the center 
(March 10) shows the decrease of the upper cutoff to ~ 470 Hz 
after three days of low activity (K < 2), and finally the one 
on the right shows the recovery of ?he upper cutoff to ~ 550 Hz 
one day after a sudden commencement. Similar decreases of the 
upper cutoff frequency have been observed during other calm 
periods (around November 18 and December 13, 1967). 

Table 1 shows the average values of upper cutoff from Rosman and 
Santiago data during the months of January and March 1968. These data 
were selected from the cases of sharpest upper cutoff, excluding the 
periods of extremely low magnetic activity, such as the days preceding 
March 10 (see Fig. 29). The table also shows the standard deviations 
(a) and solar zenith angles (^) for two approximately conjugate points 
(±35° dipole latitude). The months of January and March were chosen in 
order to provide a variety of ionospheric illumination conditions at both 
stations (as can be seen by the values of y) . The closeness between the 
resulting average values of upper cutoff frequency plus the fact that 
near Santiago the geomagnetic field intensity is ~ 25% lower than near 
Rosman indicate a lack of control of the upper cutoff frequency by the 
local ionosphere. Although the uncertainties involved in the determination 
of the cutoff frequency from the spectrograms are of the order of 3%, 
the average values for Rosman were consistently higher than those for 
Santiago (2.9% in January and 3.7% in March). This small difference 
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SNT 0G0-4 BAND-LIMITED HISS(BLH) 



FIGURE 29. DECREASE IN THE UPPER CUTOFF OF THE BLH AFTER PERIODS OF LOW MAGNETIC ACTIVITY 
AND SUBSEQUENT RECOVERY AFTER A SUDDEN COMMENCEMENT (MARCH 11). Notice the 
time lag of approximately one day between the sudden commencement and the 
recovery of the upper cutoff toward its normal value ( ~ 590 Hz) . 



may be attributed to the slight increase in upper cutoff with latitude, 
which will be discussed later (Fig. 36). 


TABLE 1. UPPER CUTOFF FREQUENCY OF THE BLH . 

January 1968 (1330 LT) March 1968 (0800 LT) 


Station 

Average Value 
(Hz) 

a 

(Hz) 

X 

No. 

Cases 

Average Value 
(Hz) 

0 

(Hz) 

X 

No . 
Cases 




IHHHI 






ROS 

582 

■ 

°o 

13 

591 

12 

63° 

n 

SNT 

565 

19 

32° 

13 

569 

15 

0 

o 

12 


As mentioned in the introduction some measurements of the amplitude 
of the signals were made using data from the sweeping receivers of OGO 4. 
These receivers cover three consecutive bands. Band 1, used here, has a 
bandwidth of 40 Hz and sweeps through 256 steps in the frequency range 
0.15 to 1.54 kHz every 73.6 seconds in the low bit rate (4 kb/sec). This 
information is recorded on magnetic tape aboard the satellite and can be 
transmitted to the ground on command from the telemetry stations, providing 
measurements along a full orbit. For more details on the experiment see 
Rorden et al. [1966]. The quantity measured is the magnetic field B of 
the waves and Fig. 30 shows an example of such a measurement for the pass 
illustrated in Fig. 26. The upper panel shows a portion of the spectro- 
gram of Fig. 26 containing, the erosion on the northern hemisphere side. 

The time scale here is expanded in relation to Fig. 26 and the diagonal 
straight line across the spectrum represents the sweep of the band 1 
sweeping receiver shown in the lower panel. The corresponding amplitude, 
in db below ly, versus frequency shows the variation in signal strength 
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QUI 0G0-4 18 FEB 1968 
1556 UT 



f (kHz) 


FIGURE 30. MEASUREMENT OF THE AMPLITUDES OF THE BLH AND THE 
BACKGROUND HISS IN THE REGION OF THE EQUATORIAL 
EROSION. The quantity being measured is the wave 
magnetic field . 
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across the two cutoffs (A and B) and the erosion (C) . The spikes 
identified as interference appear usually at multiples of 400 Hz and 
are caused by inverters aboard the satellite. Along the sweep reproduced 
in the figure, the interference at 400 Hz is not seen. However, in the 
spectrogram, the lines at 400 Hz and 800 Hz can be easily identified. 

The noise level of the band 1 sweeping receiver varies from -60 dBy 
(2.5 x 10 y /Hz) at 200 Hz to -97 dBy (5 X 10 y / Hz) at 1500 Hz. The 
small peak at 200 Hz seems to be real since it is only observed at low 
latitudes, but it has not been interpreted yet. 

In order to have a sample of how the amplitude of- the noise varies 
with latitude, measurements were made along four orbits covering a range 
of +84° to -84° invariant latitude. Two of the passes were on February 18 

and the other two on February 22, 1968. The result, plotted in units of 

2 , 

y /Hz versus invariant latitude appears in Fig. 31a, b. For each of the 
orbits, the amplitudes at three sets of frequencies were chosen for 
illustration: 600 Hz (immediately above the BLH) represented by dashed 

lines, 1000 Hz represented by continuous lines, and the frequency of 
maximum intensity in the BLH represented by dash-dot lines. Three of 
the four passes in Fig. 31 do not show recordings of amplitude for the 
whole range of latitudes covered by the satellite due to changes in the 
mode of operation of the spacecraft. 

The equatorial erosion is clearly visible and amounts to attenuations 
of 20 to 33 dB (at 600 Hz) between two consecutive sweeps, corresponding 
to a change in latitude of about 4°. The maximum intensities at 600 Hz 
and 1000 Hz occur in the vicinity of ±60° invariant latitude and may be 
of the order of 10 _4 y 2 /Hz. Beyond ±70° the intensities decrease quite 
rapidly and at approximately ±84° the background noise levels of the 
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peak between -40 o and -50 invariant latitude, but the hiss above the BLH 
peaks around ±60 invariant latitude. The gaps in the latitude coverage 
* changes in the operating mode of the satellite. 



receiver are attained. In the range of 0 to 3 kHz and up to ~ 55° 
latitude the BLH is the strongest signal observed:.. 

In Fig. 31a the maximum intensity of the BLH for: the pass of 22 

February in the southern hemisphere increases' from -30° to -46° invariant 

-4 2/ 

latitude where it peaks at ~ 2 x 10 y /Hz. A sample of the range of 

intensities at 1.0 kHz and at the frequency of maximum amplitude in the 

BLH as a function of local time was taken from 14 0G0-4 passes between 
27 December 1967 and 15 April 1968 and is shown in Fig. 32. The four 
groups of values correspond to longitude ranges of -40° to -84° (values 
around 1000 LT and 2200 LT) and 115° to'128° (values around 0400 LT and 
1600 LT) . . The four passes of Fig. 31a, b are represented in the group 
around 1000 LT in Fig. 32. The amplitudes at 1.0 kHz were taken in the 

same latitude range as the BLH in each pass. This range is limited by 

the merging of the increasing lower cutoff and the essentially steady 
upper cutoff, already illustrated in Fig. 28. For a given satellite 
altitude the lower cutoff frequency increases with latitude due to the 
corresponding increase in local magnetic field. (At higher altitudes the 
lower cutoff should occur at lower frequencies and the observation of the 
BLH would then extend to higher latitudes.) 

C. DISCUSSION 

The symmetry of both the erosion and the upper cutoff in relation 
to the magnetic equator suggests that the BLH is generated near the 
equatorial plane. Since the above properties are independent of local 
ionospheric conditions and satellite altitude, . the source is expected to 
be quite far away from the 0G0-4 orbit. 

A comparison with the observations and conclusions of previous papers 
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MAXIMUM OF BLH 



FIGURE 32. VARIATION OF THE HISS INTENSITY WITH LOCAL TIME. Both the BLH and the background 
hiss seem to peak around 1000 LT. 




on ELF hiss is restricted by the fact that most of them deal with high 
latitude and/or high altitude hiss. Guthart et al. [1968] interpreted 
the ELF hiss observed on OGO-2 in terms of electrostatic waves. According 

to their theory "...The electric field E must be coupled to the loop 

\ ..... • 

antenna, sensitive to magnetic flux density, through its cross product 
with the satellite velocity V giving the magnetic field B as follows: 

"4 — * “ 4/2 

B = V X E/c 

I T =- ■> 

where c is the velocity of light. However, for an amplitude of B 

of -50 dBy (in a 40 Hz bandwidth) the electric field required is 40 volts/m 

-4 2/ 

On the other hand, our data indicate maximum intensities of 2 x 10 y /Hz 
(-21 dBy in the same bandwidth, see Fig. 31a), which would require, in 
their theory, electric fields of the order of 1000 volts/m. Since 
electric fields of this order have not been directly observed, it seems 
appropriate to reject that interpretation. 

Gurnett and Burns [1968] report, observations of ELF hiss for in- 
variant latitudes larger than 35° and from the description of the lower 
cutoff mechanism conclude that the source of this kind of emission must 
be located above 3000 km altitude. Although these authors did not 
mention the occurrence of the BLH, examination of their Fig. 3 suggests 
its presence below ~ 600 Hz . 

Using data from 0G0 3, Russell et al. [1969] studied the spatial 
extent of ELF noise and concluded that steady noise is most common on 
the dayside of the magnetosphere within L = 6. Their strongest amplitudes 
of steady signal (> 4 volts in the 100, 300 or 800-Hz channels) showed a 
marked day-night asymmetry and were seen more frequently within L = 5, 
except between ~ 40 and 50° diople latitude, where high rates of occurrence 
were found as far out as L = 9 or 10. These authors did not report 


SEL-71-055 


- 66 



measurements below L = 2. From the Intensity of these signals they 
concluded that "on the dayside at low L values, the steady signals 
may be generated almost anywhere along a field line up to 40° magnetic 
latitude but at high L values on the dayside the signals must be 
generated in the range 40° to 50°. " 

Reporting on data from 0G0 1, Dunckel and Helliwell [1969] observed 

ELF hiss of relatively constant intensity from 0.3 kHz to as high as 

3 kHz, and from L = 2.5 up to the vicinity of the plasmapause (L = 4). 

Maximum intensities in the vicinity of the equatorial plane occur between 

06 and 16 LMT from L = 2.5 up to approximately L = 4, extending to 

L = 5.5 in the sector between 10 and 12 LMT. These maximum intensities 

_0 2 — 

are 16 dB above 10 y Hz . Based on the variation of upper cutoff with 
latitude, these authors concluded that the signals are generated close to 
the equatorial plane. 

Studying the direction of the poynting flux of ELF hiss with the 
Injun 5 satellite in the region 677 to 2528 km altitude and 35° to 75° 
invariant latitude (in the northern hemisphere), Mosier [1971] concluded 
that the hiss is generated outside the plasmasphere and propagates across 
the plasmspause at altitudes probably below 2500 km. This conclusion is 
based on his remark that " . . .downgoing ELF hiss was observed over the 
entire region of altitude-invariant latitude space under study, whereas 
the highest latitude at which the net poynting flux of ELF hiss was 
directed upward is approximately 60° invariant latitude." 

The possibilities offered by an equatorial source can be tested 
using raytracing techniques. The results, as described in the following 
sections, are used to interpret the several characteristics of the BLH, 
the observations of ELF hiss outside the plasmapause between 40 and 50 
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reported by Russell et 1. [1969] and also some, of thepoynting flux 
measurements of Mosier [1971]. In these ray tracings, the ionosphere 
model used is in diffusive. equilibrium along magnetic field lines and 
describes quite well the daytime equatorial anomaly in the electron 
density. It is the same model used previously in Chapters 2 and 3. 

1 . . Accessibility to Low Latitudes 

Figures 33a, b,c show raypaths for a frequency of 600 Hz starting 
at L = 3 in the equatorial plane, for several values of the initial 
wave normal angle The wave normal angle lj; is measured from the 

magnetic field direction so that |t|i | < 90°, and is positive when the 
wave normal points toward lower L shells. The field line represnted at 
L = 4 marks the position of a simulated plasmapause obtained by smoothly 
decreasing the electron density by a factor of 10 over a distance corre- 
sponding to AL = 0.23. 

The lower hybrid resonance (LHR) frequency at L = 3 in the equatorial 
plane is ~ 750 Hz and so the refractive index surface at 600 Hz is closed. 
After the waves leave the equatorial plane, the behavior of the wave 
normal angle depends on its initial value, in a manner consistent with 
that shown by Thorne and Kennel [1967]. In the present case if \|i < 50° 

(Fig. 33a), \|j decreases continuously until it reaches -90°, at which 

point the rays are refracted upward. 

For \j/ . > 55°, ij) first decreases, reaches a minimum value and 
then increases again. In addition, with increasing the minimum 

value of ijf also increases, and the point at which this minimum occurs 
approaches the equator. 

If \Jj . > 70° the minimum wave normal angle is reached so close 
to the equator that tj/ may be able to increase up to +90 , when the rays 
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FIGURE 33a. RAY PATHS OF 600 Hz WAVES STARTING IN THE 
EQUATORIAL PLANE AT L = 3 FOR DIFFERENT 
VALUES OF THE INITIAL WAVE NORMAL ANGLE ty . 
The field line at L = 4 marks the position of a simulated 
plasmapause . The rays corresponding to tji ^ between 40° 
and 50° escape from the plasmasphere . 
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FIGURE 33b, ADDITIONAL RAYS FOR .A. v BETWEEN 70° AND 
o Y 1 

85 . The arrows indicate wave normals at 

several points along the rays corresponding 

to i|j equal to 70 and 75°. 


FIGURE 33c. RAYS CORRESPONDING TO i]> . EQUAL TO 86° AND 
89°. In order to reach latitudes below 30° 
at OGO-4 altitudes, the initial wave normal 
angles have to be larger than 75°. 
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will be refracted inward producing the kind of raypaths illustrated in 
Fig. 33b. For these raypaths, each time the equatorial plane is crossed 
the wave normal moves closer to the geomagnetic field direction (ty decreases) 
and it may happen that after a few zig-zags, the equatorial value of \jj 
will be so low as to cause the next refraction to be upward (as for . = 70° 

in Fig. 33b after the first refraction). If the initial wave normal 
angle is sufficiently large (ty > 75° in the figure) the rays will bounce 
back and forth across the equatorial plane on their way down and these 
are the only rays capable of reaching latitudes less than ~ 30° (Figure 
33b, c; the wave normals corresponding to i|f = 70° and 75° are indicated 
by arrows in Fig. 33b). Although the initial wave normal angles have to 
be large, the wave normal angles at low altitudes will be influenced by 
the ionospheric density gradients and can assume different values. 

2 . Erosion 

In order to understand the distribution of energy around the 
equator, ray tracings were made at intervals Ai|/ ^ = 0.1° from 85.5° up 
to 89°. Two frequencies were used, 500 Hz and 600 Hz and the results are 
plotted in Fig. 34 as latitudes of arrival (at 700 km altitude) vs initial 
wave normal angles. The points distribute themselves in two groups on 
both sides of the equator and for a uniform increase in \J;^ there are 
sudden jumps from one hemisphere to the other. This defocusing effect 
produces an attenuation of signals across the equator, as observed in 
the satellite recordings. It. is possible to make a rough estimate of 
the relative intensity of the signals by assuming that the intensity is 
inversely proportional to the spacing between adjacent points. For the 
distribution points in Fig. 34 the result is shown in Fig. 35, where 
intensities (in an arbitrary linear scale) are plotted vs the absolute 
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RELATIVE INTENSITY 


FIGURE 35. RELATIVE INTENSITIES OBTAINED FROM THE DISTRIBUTION 
OF POINTS OF FIG. 34 VS DIPOLE LATITUDE. The 
variation in intensity is consistent with the 
equatorial erosion. 
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value of dipole latitude. The irregularities in the resulting curves are 
due to the utilization of the discrete (and small) number of points from 
Fig. 34. If a larger number of points had been used, the intensity vs 
latitude curves of Fig. 35 would be smoother and also their slope would 
be steeper. Nonetheless, from ~ 10° to ~ 17° the intensity for 500 Hz 
is consistently higher than the one for 600 Hz, and a receiver with a 
constant threshold sensitivity in this range of frequencies (as in 0G0 4) 
would observe an erosion as in Fig. 26. 

3 . Upper Cutoff 

Another characteristic of the BLH reported here is the apparent 
constancy of the upper cutoff frequency observed in spectrograms from 
the broadband receivers. However, when examining the amplitude data from 
the sweeping receivers one sometimes notices a slight displacement (vari- 
able from pass to pass) of the upper cutoff toward higher frequencies 
with increasing latitudes. This displacement is, however, much smaller 
than the corresponding variation of the lower cutoff. To illustrate 
that displacement the outputs of the band 1 sweeping receiver corresponding 
to parts of three OGO-4 passes are shown in Fig. 36. For each of the 
passes the consecutive sweeps are shown closely spaced in order to enhance 
their differences. Latitudes increase toward the bottom of the figure 
and the mentioned displacement is quite noticeable, mainly in the pass 
of 22 February 2200 LT. The sharpness of the upper cutoff is also 
variable and for the pass of 27 December it has a slope ranging from 8 
to 15 dB,/l00 Hz. In the sample of 14 passes examined, this slope varied 
between 8 and 20 dB/l00 Hz, while the lower cutoff was usually steeper 
(between 15 and 50 dB/l00 Hz) . 
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0.2 0.6 l'.0 l'.4 kHz 

27 DEC '67 
1536 LT 


1000 LT 




22 FEB '68 
2200 LT 


FIGURE 36. REPRODUCTION OF BAND 1 SWEEPING RECEIVER OUTPUTS 

FOR PART OF THREE OGO-4 PASSES. Subsequent sweeps 
are reproduced closely spaced in order to enhance 
the variations in shape of the BLH. Latitudes increase to- 
ward the bottom of the page and one can notice a slight 
displacement of the upper cutoff toward higher frequencies 
with increasing latitudes. 
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The upper cutoff of the BLH can also be understood if the source 


is located in the equatorial plane. The ray tracings of Fig. 33, starting 
at L =3, were intended to illustrate the general behavior of rays in- 
side the plasmasphere and they indicate the need for large initial wave 
normal angles in order to reach latitudes less than ~ 30° at low 
altitudes. If the source altitude is varied, however, the behavior of 
rays will not differ much from those illustrated. In the discussion of 
the erosion effect it was mentioned that the refractive index surfaces 
corresponding to the waves propagating toward lower L shells in Fig. 33 
were closed. (The LHR frequency at L = 3 in the equatorial plane is 
~ 750 Hz.) It is necessary that the surfaces be closed in order to produce 
the sudden refraction characteristic of the zig-zag ray paths, otherwise 
the wave normal angle could never attain 90°. When the frequency is 
increased past the equatorial or initial value of the LHR, the refractive 
index surfaces become open in that region but it is still possible to 
generate zig-zag ray paths provided the ratio is not much 

greater than unity. This is so because the rays will be moving away from 
the quatorial plane approximately along magnetic field lines and the 
local magnetic field intensity will increase quite rapidly (and the f 

LHR 

with it). Soon a point will be reached where the local f becomes 

LHR 

larger than the wave frequency and the refractive index surfaces will be 
closed again. From the results of ray tracings with frequencies that 
obey the above condition it was found that when the frequency is increased: 

1. There is an increase in the minimum value of initial wave 
normal angle (ty.) necessary to generate a zig-zag raypath 
down to 0G0-4 altitudes. 

2. For this minimum value of tjt . , the latitude of arrival of 
the rays also increases slightly. 


77 


SEL-71-055 



However, as the frequency is increased above the LHR frequency the 

resonance angle (where the refractive index p. -> oa) decreases. Actually 

when ill . becomes larger than ill (defined as the value of \1 i for which 
i m 

the product pcos\)i is minimum) the rays begin, to move toward larger L 
shells and are not able to reach low altitudes and low latitudes. If 
the source of the hiss is in the equatorial plane, therefore, a low 
altitude satellite like 0G0 4 will observe an upper cutoff (somewhat 
above the LHR frequency at the generating region) that should increase 
slightly with latitude. In the examples of Fig. 33 with the origin at 
L = 3 this upper cutoff happens above 1 kHz. In order to obtain a cut- 
off close to 600 Hz the source has to be moved to a higher altitude. If 
the source is moved to L = 4 the equatorial value of the LHR frequency 
is reduced! to approximately 320 Hz. Some rays were actually traced for 
a frequency of 600 Hz starting at that altitude and it was found impossible 
to make any rays arrive at latitudes lower than 50° and low altitudes. 

A source close to L = 4 would, therefore produce a hiss band with 
characteristics similar to the observed BLH and its upper cutoff would be 
essentially independent of the magnetic latitude and local conditions 
at the satellite. 

The observed decrease of the upper cutoff of the BLH shown in 
Fig. 29 would also be understood, in the context of the present interpre- 
tation, as an outward drift of the generation region during periods of 
very low magnetic activity. 

4 . Leakage from Plasmasphere 

Transmission of rays to regions outside the plasmasphere can occur 

‘ o o 

for initial wave normal angles i|i varying from 40 to 50 as shown in 
Fig. 33a. This result may be used to explain the observations of 
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Russell et al. [1969] of daytime occurrences of "steady noise" outside 
the plasmasphere between 40° and 50° dipole latitude. The rays with 
< 40° are refracted back to the plasmasphere by the steep density 
gradient that constitutes the plasmapause and are trapped inside. The 
ray path corresponding to ij; = 55° is bent downward, also by the gradient 
at the plasmapause, and if ^ > 55° the rays are confined within the 

plasmasphere. Even when the rays were started at higher altitudes, as 
mentioned in the previous section, the leakage to the outside across 
the plasmapause could still be observed around ~ 45° dipole latitude. 

5 . Direction of Poynting Flux 

In his study of the direction of poynting flux using data from the 
Injun 5 satellite, Mosier [1971] reported the occurrence of downgoing 
ELF hiss over a range of invariant latitudes from ~ 40° to more than 
70° (across the plasmapause) . The highest invariant latitude at which 
he observed the poynting flux of the hiss directed upward was about 60° 
(approximately the plasmapause location) and during daytime these 
occurrences were always above 1500 km. Again, from the ray paths in 

Fig. 33a, it can be concluded that below ~ 3000 km and outside the 

plasmasphere the rays are always propagating downward, but in the plasma- 
sphere rays propagating both downward and upward are encountered (Fig. 

33b). If 70° < i|r < 75° for example, some rays were traced and 

refracted upward at altitudes as low as 2500 km. With a nighttime 

density model this altitude was easily lowered to 1200 km, using fre- 
quencies as high as 900 Hz. As mentioned before, the wave normal angles 
at low altitudes are sensitive to density gradients and are, therefore, 
dependent on the ,model used . 
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D. GENERATION MECHANISMS 

The characteristics of the BLH presented and discussed in previous 
sections suggest a source location near the equatorial plane in the 
vicinity of L = 4 and the necessity of large initial wave normal angles. 
In this section we discuss briefly two possible generation mechanisms 
and make a simple estimate of the energy levels and particle numbers that 
can reproduce the measured values of hiss intensity. 

The two most popular radiation mechanisms to explain VLF and ELF 
emissions from the magnetosphere are Cerenkov and Doppler-shifted 
cyclotron (DSC) radiation from electrons trapped in the earth’s magnetic 
field (Ellis [1957]; Gallet [1959]; Gallet and Helliwell [1959]; Dowden 
[1962,1963]; Mackenzie [1963,1967]; Bell and Buneman [1964]; Brice [1964]; 
Liemohn [1965]; Kennel [1966]; Kennel and Thorne [1967]; Helliwell [1967]; 
Mansfield [1967]; Jorgensen [1968]; Melrose [1968]; Russell and Thorne 
[1970]; Trulsen and Fejer [1970]). Without going into details of the 
radiation processes (see above references) the conditions for Cerenkov 
(or Landau) and DSC resonance between waves and electrons can be described 
by the following expressions [Liemohn, 1965]: 

pcostk = — for Cerenkov (or Landau) resonance 
Pii 

1 f H 

ucosili = — - - (1 - m — ) for DSC resonance 

* Pn Y f 

In these expressions p is the wave refractive index, iji is the wave 

normal angle of the emitted radiation, j3 () is the parallel component of 

the normalized particle velocity (|3 = — ), f„ is the electron gyro- 

C 2 - 1/2 

frequency, f is the frequency of the radiation, y = (1 - (3 ) 
accounts for relativistic effects and m is an integer (the case m = 0 
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H 

corresponds to the Cerenkov emission). When — - > 1, the emission is 

yf 

either in the forward or the backward hemisphere depending on the sign 
of the integer m. If m < 0 the emission is in the forward hemisphere 
(c stjj and g of the same sign) and if m > 0 the emission is in the 
backward hemisphere (costy and (3 of opposite signs) as illustrated in 
Fig. 37. 



FIGURE 37. RELATIVE POSITIONS OF THE PARALLEL COMPONENT 
OF THE VELOCITY OF THE EMITTING ELECTRON AND 
THE WAVE NORMAL, IN A MAGNETIC MERIDIAN PLAN, 

FOR THE CASE OF CERENKOV (OR LANDAU) RESONANCE (m = 0) 

AND DOPPLER SHIFTED CYCLOTRON RESONANCES (m > 0) . 

For the range of frequencies encountered in the BLH and the source 
location discussed here (near L = 4) the product pcostji can be as low 
as 18 due to the large wave normal angles involved. If one considers the 
Cerenkov process the above expression- for the condition of resonance 
indicates electron energies of the order of 1 kev On the other hand 
the DSC mechanism will involve electron energies near 1 Mev. 

Expressions for computing the power emitted by a single electron 
on the Cerenkov and DSC mechanisms can be found elsewhere (see Mansfield 
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[1967]; Melrose [1968]; Trulsen and Fe jer [1970]) . The quantity computed 
is the power radiated per frequency interval (dP/df) as a function of the 
particle velocity components v and v respectively parallel and 
perpendicular to the earth's magnetic field, and the characteristics of 
the medium. For a frequency of 400 Hz and a wave normal <|j • = 87° in 
the equatorial plane at L = 3.8 we found: 

1. For the Cerenkov mechanism: 

electron energy: 0.56 kev 

^ i 4. x 10 watts/Hz 

2. For the DSC mechanism: 

electron energy: 780 kev 

dP -34 / 

— = 1. X 10 watts/Hz 

The emitted power per electron is of the same order of magnitude in both 

mechanism. Frank [1968] indicates omni-directional electron fluxes 

9 / 2 

(0.3 to 3 kev) of 2 x 10 el/ cm sec around L = 4. For the high energy 

particles (E > 500 kev) Vernov et al. [1969] indicate intensities of 

about 6 x 10 el/cm sec. Comparing these two intensities and considering 

the differences in particle velocities, the number density of the low 

4 4 

energy electrons will be from .10 to 3 X 10 times larger than that of 
the high energy electrons. 

Taking the above value of dP/df obtained for the Cerenkov emission 
mechanism, the power flux emitted by a single electron can be computed. 

By comparing this result with the power fluxes measured by 0G0 4, an 
estimation of the order of magnitude of the number of emitting electrons 
can be made. In a frequency interval Af each electron radiates through 
an annulus subtending a wave normal angle Aty around \ji but only a small 
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segment of this annulus, pointing toward the earth, will contribute to 
the power received by the satellite. Taking for instance a segment of 
2° (that is, considering wave normals slightly off the meridian plane) 
and making the conservative assumption that only rays starting inside 
this segment will propagate down to OGO-4 altitudes it is possible to 
estimate the area AA illuminated by this segment of spherical zone at 
OGO-4 altitudes for a given frequency band Af and the corresponding 
Ail) • (The rays off the meridian plane were traced utilizing a three- 
dimensional computer ray tracing program developed by F. Walter.) 

Taking for Af the bandwidth of the band 1 sweeping receiver (40 Hz) 
the area AA is found to be 

11 2 

AA A 3 .6 X 10 m 

which corresponds to an emitted power flux per electron (S ) of 

La 

— 46 2 

S = 2.5 X 10 watts/m 

E 

—6 2 

From the intensity found in the BLH at low latitudes ( ~ 5 X 10 y /Hz, 
Fig. 31a, b) and assuming a value of 50 for the refractive index, the 
measured power flux on the same 40-Hz bandwidth is found to be of the 
order of : 

-10 / 2 
S = 4 x 10 watts/m 
M 

Assuming that the hiss is generated incoherently, the ratio of these 
power fluxes gives the number of emitting electrons: 

S 

N = — £ 1.6 X 10 electrons 
S E 

It is easily verified that this is an excessively large number, since 
even for a density of 10 el/cm , it corresponds to a volume greater 
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than ten thousand magnetospheres. Therefore, some form of coherency in 
the emission generation, or wave growth, is needed to produce the observed 
power fluxes. It. can easily be verified that similar requirements also 
hold for the Doppler-shifted cyclotron generation.. However, a detailed 
study of the possibilities of wave growth. in the DSC mechanism is outside 
the scope of this work. An excellent treatment of the problem can be 
found in Russell and Thorne [1970], although these authors give emphasis 
to small wave normal angles and frequencies higher than in the BLH. 

A greater knowledge of the ;low energy electron fluxes in the vicinities 
of the plasmapause and their drifts with magnetic activity could provide 
a test for comparison of the two emission mechanisms. 
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V. CONCLUSION 


A. INFLUENCE OF THE MEDIUM ON d\ji/ds 

The study made in Chapter 2 of the problem of orientation of the 
wave normals established the relative importance of several factors such 
as the electron and ion densities, the magnitude and inclination of the 
earth's magnetic field on the rate of variation of the wave normal angle 
ij . In the examples given, the frequency of the waves and the range of 
altitudes chosen put emphasis on the effects of the hydrogen ions, 
which are the most important in connection with the ion whistlers 
studied in Chapter 3. 

The inclination of the magnetic field (dipole) plays a small role 
in the altitude range examined (450 to 2000 km) and its effect can 
safely be neglected, when compared to the other two terms of di|//ds. 

The effects of these two terms are mode dependent and can be summarized 
as follows: 

1 . The Particle Density Term 

For the slow mode, it is possible to separate two altitude ranges 

where different contributions dominate. Below this transition altitude 

the important contribution is that of the electrons and the wave normal 

will turn toward the direction of VN . The transition altitude is 

e 

situated close to the minimum of p(r). It approaches the crossover 

altitude from below as the wave normal angle tends to zero. Above the 

transition, the important contribution is that of the hydrogen ions and 

the wave normal will turn toward the direction of VN +• For any wave 

n 

normal angle different from zero, the particle density term goes to 
infinity at resonance. In Appendix A, curves of dty/ds are plotted for 
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the same conditions as in Figs. 8 to 11, but for a wave normal angle of 
89°. The decrease in the resonance altitude toward the two-ion resonance 
can be seen in those curves . 


For the fast mode, which exists only above the two-ion cutoff 

altitude, there are also two regions where different contributions 

dominate. The transition altitude approaches the crossover altitude 

from above when the wave normal angle decreases to zero. Below the 

transition the dominant contribution comes from the hydrogen ions and 

the wave normal turns toward the direction of VN^. Above that altitude 

the principal contribution is from the electrons and the wave normal 

turns toward the direction of VN . 

e 

■4 

2 . The Term Due to the Magnitude of B 

For the slow mode and small wave normal angles, this term becomes 
important between the crossover and the proton gyrofrequency altitudes, 
where the refractive index varies appreciably with B. Its magnitude 
grows to infinity at resonance, but unlike the particle density term 
the sign is conserved. 


In the fast mode this term is important either near the two-ion 

cutoff altitude (where — -* ®) or at large altitudes where the particle 

u OB . 


density term becomes negligible. 


B. ELF HISS .. . 

Analysis of ELF and VLF data from 0G0 4 from September 1967 to 
April 1968 revealed some features of the ELF hiss that has not been 
previously identified (until recently, the only special feature of the 
hiss that had been identified was the lower cutoff). In particular, 
the hiss band below approximately 600 Hz (called BLH, for band limited 


ELF hiss) presents the following features: 
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1. A lower cutoff frequency increasing with latitude, typical 
of downward propagating hiss and attributed to the two-ion 
cutoff in the dispersion relation [Gurnett and Burns, 1968]. 

2. An upper cutoff frequency close to 600 Hz which seems to be 
essentially indepeneent of latitude above ~ 15° or 20° (cf. 
Fig. 26) (sometimes it can show a slight increase toward higher 
latitudes — cf. Fig. 36). After a prolonged period of low 
magnetic activity (K p < 2) the upper cutoff frequency may 
decrease to approximately 420 Hz but recovers about one day 
after a sharp increase in K p or a sudden commencement 

(Fig. 29). (The data examined did not include any period of 
intense magnetic activity, K p > 6.) 

3. The upper cutoff decreases in frequency toward the equator 
near 10° dipole latitude (equatorial erosion) (Fig . 26). The 
attenuation rates in the region of erosion are in the range 
of 20 to 33 db at 600 Hz for a latitude change of ~ 4° 

(Fig. 31a, b). 

4. The latitude of occurrence extends from ~ 10° up to ~ 55° 
dipole latitude, where the essentially constant upper cutoff 
and the increasing lower cutoff merge (cf. Fig. 28). This 
limiting latitude depends on the local magnetic field intensity 
and the satellite altitude. 

5. The most frequent times of occurrence lie between 0600 and 
2200 LT (when it showed up in almost every OGO-4 pass), 
although some occurrences have been found between 0200 and 
0500 LT (with less intensity). 

6. The local ionosphere has no apparent affect upon the upper 
cutoff frequency, which is virtually the same at conjugate 
points . 

Between the latitudes of ~ 10° and ~ 55° the BLH contains the 

strongest signals in the frequency range of 0 to 3 kHz. Its intensity 

peaks between 40° and 50° dipole latitude and may reach 2 X 10 ^ y^/Hz 

during daytime ( ~ 1000 LT, cf . Fig. 31b). The decrease and smearing of 

o 

the lower cutoff of the background hiss, observed above 58 invariant 
latitude, is an indication of the light ion trough. 

The characteristics of the BLH suggest a source location in the 
equatorial region near L = 4 and the need for large initial wave normal 
angles. This source location provides an explanation of: 

1. The presence of the BLH at low latitudes and the equatorial 
erosion; 
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2. The nearly constant upper cutoff of BLH; 

3. ' The daytime observation of hiss outside the plasmapause between 

40 and 50 latitude reported by Russell et al . [1969]; 

4. Some observations of poynting flux direction reported by 
Mosier [1971] . 

The Cerenkov and Doppler-shifted cyclotron (DSC) generation 
mechanisms from trapped electrons are briefly examined as possible sources 
of the BLH. For the frequency range in question the electron resonance 
energies are of the order of 1 kev for the Cerenkov emission and close 
to 1 Mev for the DSC emission and the power radiated by a single electron 
is comparable in both mechanisms. In either case, however, some degree 
of coherency or amplification is necessary in order to explain the 
observed power fluxes. At this point it is not clear which one of the 
two emission processes is responsible for the BLH. 

The interpretation of the upper cutoff is consistent with a decrease 
in the population of the generating particles inside L = 3.5 or L = 4 
and an outward movement of this generating region at times of very low 
magnetic activity. Therefore, a correlation study between the variations 
in the upper cutoff of the BLH and movements of energetic electrons 
could aid in finding which particles are responsible for the BLH. 

The conclusions of the present work do not preclude the possibility 
that the same mechanism that is responsible for the BLH may also be 
generating the background hiss. As a matter of fact, the equatorial 
erosion seen in the background hiss in Fig. 26 suggests exactly this. 

In this context, differences in particle fluxes at different altitudes 
could be the important factor in explaining the different hiss intensities 
above and below ~ 600 Hz . 
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C. SUGGESTIONS FOR FUTURE WORK 

Based on the analysis and interpretation of OGO-4 ELF hiss obser- 
vations (and the BLH in particular) presented here, the following topics 
are recommended for further study: 

1. The variation of the upper cutoff of the BLH with magnetic 
activity. The study made in Chapter 4 was limited by the 
availability of good data at times of very strong and very 
weak activity. The question of the time lag between particular 
magnetic fluctuations and the corresponding effects on the 
upper cutoff needs to be ascertained more accurately. 

2. The variations of the amplitude of the BLH with local time. 

The results shown in Fig. 32 are preliminary (only 14 satellite 
passes were examined). Also the variations of amplitude of 

the BLH with magnetic activity should be examined. 

3. Correlation of the properties of ELF hiss with fluxes and move- 
ments of energetic particles (electrons) near the plasmapause 
This study can help clarify the question of which particles 
are responsible for the observed hiss bands . 

A great deal can be learned about the above problems by an extensive 
analysis, on a global scale, of amplitude data recorded by the PCM 
receivers of OGO 4. 

The possibilities of physical bunching or emission coherency 
should also be investigated together with an examination of the required 
lifetime of the emitting particles. 

4. Correlation of the smearing and disappearance of the lower 
cutoff of the ELF hiss at higher latitudes and the plasmapause 
location. If a good correlation can be found the above 
phenomenon could provide a very easy and efficient way of 
determining the plasmapause location. 
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APPENDIX A. DERIVATION OF AN EXPRESSION FOR difr/ds 


In what follows the medium is assumed to be lossless, so that the 
energy flow will be along the group velocity vector [Stix, 1962]. This 
direction is also called ray direction and will be denoted by s. The 
velocity of the wave fronts in the ray direction, or ray velocity has a 
magnitude related to the magnitude of the phase velocity by the expression: 


~ ds 
V R ~ d t 


v _ c 
cosa ~ pcosa 


(A . 1 ) 


where a is the angle between the direction of the phase and ray 
velocities (see Fig. 38), and can be expressed by: 


tana 


1 

^ ^ r , 6 


(A. 2) 


In the present work the earth's magnetic field (see Stix [1962]) is 

assumed to be a dipole. A further simplification will be made by only 

dealing with waves propagating in magnetic meridian planes. Polar 

coordinates are therefore very convenient, and the polar axis is chosen 

along the axis of the dipole. So at the point P(r,9), B is the 

magnetic field vector and v and v are the phase and ray velocity 

R 

vectors, respectively. The components of the ray velocity are: 

^ = v cos(6 + a) (A. 3) 

d t R 

r cTE" = V R sin ^ + (A. 4) 

V 

— = — cosa(sin6 + tana cos5) 
d t r 

and by virtue of Eqs . (A.l) and (A. 2) 


d0 

dt 



p sin6 - ( - 

It) cos6 

\ 

■'r,Q J 
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(A. 5) 



/ 

POLAR AXIS 



FIGURE 38. POLAR COORDINATE SYSTEM. All angles are measured 
clockwise from the references shown. 


In Fig. 39, AB represents a small extent of a wave front at 
instant t, and CD is part of the same wave front after an interval 
At. If the phase speed of point A is v, after the interval At 
this point will be displaced by an amount: 


AC = v dt 


Point B will be displaced by: 


— d v 

BD = V + rr d£ 
d-t ik 


dt 
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POLAR AXIS 



FIGURE 39. GEOMETRY OF A PROPAGATING WAVE FRONT. After 

a time interval At the wave front AB 
moves to CD . 


Therefore: 

tan(dlj)) = d<|) = - (^) dt 

Q 

and using the refractive index p = — : 

^ = i) (A. 6) 

dt 2 \d IIl 

u 9 

It is important to note that the derivative of p in Eq . (A. 6) is 
obtained by keeping the wave normal direction fixed in space. 

Because p = p [r, 0 ,ijj (r , 0 , 9) ] , the term (dp/dt)^ may be written 

as : 
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(A. 7) 


but (Jf = <() - g 


and for a dipole field 


a*\ 
3e JL 




= “\3eJ 

\ /r 


(i) 


e,<j) 


= o 


From Fig. 39 it can also be seen that: 


dr 

dt 


— = - sing 


and 


de 

61 


1 . 
— COSQ 
r 


Therefore, Eq. (A. 7) can be written as: 


feL = ~ sin6 (If) + 7 cos6 [(le) " (^) (If) 

x v '94 L\ 'r 4 \ T 'r,e \ 'r J 


Now if the unit vectors k and k are defined as: 

1 


k = r cos6 + 9 sin6 


k = r sin6 + 0 cosfi 


k and k^ are, respectively along the perpendicular to the wave 
direction. Consequently, Eq. (A. 11) becomes: 

(«l = <% V ■ i c ° s6 (ll, e (^l 


Now the rate of change of the wave normal angle will be: 

dj, _'d6 d£ 

dt “ dt dt 


and using Eq. (A.6): 


dt 


c /<V\ 

* „ 2 W* ' 


dp 

dt. 


(A. 8) 


(A. 9) 


(A. 10) 


(A. 11) 


norma 1 


(A. 12) 


(A. 13) 
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For a dipole field 


dg. /fe\ de 

dt \ ae / dt 


and combining with Eq . (A. 5): 


dg c 

dt ~ 2 

rp . 


(Dj*' -($) 


cos 5 


(A. 14) 


So, from Eqs . (A.12), (A. 13) and (A. 14), the rate of change of the wave 
normal angle can be expressed as: 


djj_ _ c 
dt ~ p 


. £ . i (i 

*§) sin6 

P X r \c 

»/ r J 


(A. 15) 


Now from Eq . (A.l): 


d^ _ d£ c_ 


dt ds p cosa 
and again for a dipole field: 


r ( ) Sin6 = 
\ / r 


Therefore, the rate of change of the wave normal angle along the ray path 
will be: 


dij(_ 


ds 


cosa 


(Vp), 


• k - Vg *k 


The gradient of p at a constant i|t may be written in terms of the 
electron and ion densities and the magnitude of the magnetic field, so: 


31 _ 


ds 


= cosa 


Z-4^-VN.-k +-*VB-k - 96-k 

L j „ 3T j x m % i. 


(A. 16) 


It is interesting to note that, if the direction of had been 

chosen to be \|j - n/2 instead of $ + n/2, the vector k^ would be 
written as: 


k = r sin6 - 0 cos6 
X 
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a na 


the 




6Ss io« 


tof 


as 


in 


Ce ^isi 


d $/c is 

St r Woul “ be; 

ds ' ~ cos « jH 

/jU Vv* *i«k 


er 


^ l 9 6 ?J < 


i + 
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APPENDIX B. MAGNETOSPHERE DENSITY MODEL 


The magnetosphere density model is composed of electrons and the 

ions of hydrogen, helium and oxygen. Electrons and ions are in diffusive 

equilibrium along lines of the earth's magnetic field (see Angerami and 

Thomas [1964]). The temperature is assumed to be the same for all 

particles. The electron density N g at a given point P(r,6) (see 

Fig. 40) depends only on its value at a reference level r^ along the 

same magnetic field line. Because of this constraint, the' latitudinal 

variations of density may be described at this reference level by a 

function n(9 ) of the magnetic colatitude 9 of the base of the field 
E E 

line in question. In other words: 

N (r ,0 ) = N • n(9 ) (B.l) 

e o o e E t 

o 

in which N is a constant, 
e 

o 

Therefore, given a point P(r,9), the corresponding colatitude 
at the base of the field line is immediately determined: 

9 - 9_(r,0) (B.2) 

EE 

and the electron density at the point will be: 

N (r, 9) = N • n(0 ) • n_(r,0 ) (B.3) 

e e e u e 

o 

where n (r,0 ) is a function describing the radial dependence (diffusive 
D E 

equilibrium) and also eventual latitudinal variations in composition and 
temperature. If the composition and temperature are to be constant with 
latitude this last function will simply be: 

V r ’V- V r> (B,4) 
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where : 


z = r (1 - — ) is the geopotential height (neglecting the 
centrifugal force) 


a = a (8„) 

o . o . E 
i 1 


h. = JES 

i m.g 


1 o 


K 

T 

g. 


"t h 

is the relative concentration of the i ion 
species at the reference level and may be variable 
with latitude 

is the scale height of the i neutral species 


is the Boltzmann constant 

is the electron and ion temperature 

is the acceleration of gravity ft r 


Due to the electrical neutrality of the medium: 


I a =1 
o . 

l l 


The ion densities will then be: 


where : 



N 

e 


a i 


-z/H. 

V e 


2 

n„ 


(B.7) 


(B.8 ) 


and also: 



Therefore, from Eq.(B.3) the derivatives of the electron density may be 


written as : 


dN 

e 

1 

Sn D 

1 

dn 

S9 e 


‘ n D 

+ 

n 


Sr“ 

3n 

e 

1 


1 

dn 

w 

CD 

CO 

w = 

n D 

3T + 

n 


ST 


(B. 9) 


(B.10) 


and for the ions: 
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(B.ll) 


1 

c)N. 

1 

, Sn 

1 e 

1 _ ^i 

N. 

= 

N cSr~ + 

a. <^r 

1 


e 

i 

1 

•H 

, <)N 

1 e 

1 _ ^i 

N. 

1 

w = 

r 5r + 

e 



( B. 12) 


The derivatives of -0 in Eqs . (B.9) and (B.10) can be obtained from 

‘ ' • * .fci •• ‘ 


Eq . (B.5) : 


3F - 


2r • (L-l) 


172 


(B.13) 


* 9 e 

sr 


tan9 • (L-l) 


1/2 


(B. 14) 


(It is important to note that 9 < */2.) 

E 

The derivatives of the function n^ appearing in Eqs. (B.9) and 
(B.10) can be obtained from Eq . (B.6): 




a. - ckT" a <39 3T" 

1 O . 

1 


(B.15) 


(B. 16) 


(B.17) 


( B. 18 ) 


In order to completely determine the model it is necessary to 

specify the functions n(0 ) and a (9_), which describe, respectively, 

E o. E 

1 

the variations of electron density and composition with latitude. Through 
out the present work, the composition is made constant with latitude, 
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e.g., a (9,,) = constant. The function that represents the latitudinal 

O . E 

1 

variations in electron density, on the other hand, is made out of the 
product of two other functions n_ and n that reproduce the so-called 
equatorial anomaly and the "knee" or plasmapause. These functions are, 


respectively : 


n = 1 + A* exp 
G 


_ 1 ^ 9 peak " 9 e | 


(B . 19) 


n K = B + C 


• erf ( J T77 L ) 


CB.20) 


and they are reproduced in Fig. 41a, b, with the values of the constants 


used in the examples. The function n is a gaussian curve, and the 

G 

function n is made with the so-called "error function." Their 
K 

derivatives which appear in Eq. (B.9) and Eq. (B.10) can be written: 

dn 
^E 

v2 ~i 


1_ 

n. 


G 

3e~ 




i_ 

n 


dn 


K 


K 




exp 


E 9 K 

2/V 


) 


(B.21) 


(B.22) 
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FIGURE 41a, b. FUNCTIONS USED TO SIMULATE THE EQUATORIAL 
ANOMALY (a) AND THE PLASMAPAUSE (b) IN THE 
MODEL IONOSPHERE. 
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APPENDIX C. DERIVATION OF EXPRESSIONS FOR dp/dN AND dp/dB 


The refractive index p satisfies the following equation: 

„ -1 


tan 


o _1 2 “I 

( 1-R/m. ) + (1-L/u ) -2 


+ 2 


(i-p/n ) -1 


= 0 

(C.l) 


where the quantities R, L and P have been defined in Chapter 2. In 
other words : 

|i = |i(R,L,P,i|i ) (C . 2 ) 


Therefore : 


dp dp dR dp dL dp dp 

diT = dii diTT + dl dN~ + Bp diT 

j j j i 


where : 


dp 


dR *j . dL S _ 

dir n . ’ . , 3n7 n . ’ d>r n . 

J J 3 3 3 3 


Consequently : 


1 ^ 1 ( 1 ^ R , 1 ^ T 4- 1 C*Li. Y \ 

4 5T = -5~U 5'jV 5 L i* [T3p X jj 


(C.3) 


(C .4 ) 


(C.5) 


Similar ly : 


dp _ dp dR dp dL 
dB _ dR c>B + dL di 


(C .6) 


where : 


dR 

35 


dR . 
v J 

E 3b~ 

J 


1 S j Y j 2 
J J 


dL 

dB 


dL . 
v J 

E ds~ 

J 


, e .Y . „ 

1 v 3 3 t 2 

' B S T~ L j 
J J 


(C .7) 


Consequently : 


1 dp _ 1 £ j Y j 
dB '* B . X . 




1 dp 2 _ 1 dp 2 \ 
p dR j p dL j / 


(C .8 ) 


3 3 

The derivatives of p with respect to R, L and P can be 
easily obtained from Eq . (C.l) above: 
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Analogously : 




(C .9) 


(C.10) 


(C.ll) 


Limits of the above derivatives at resonance (p. -« <”) . The condition 

for resonance can be found from Eq. (C.l): 

, 2 . P 

tan = - - 

where S = — (R + L) (see Stix [1962]). Therefore, at resonance (and 
for / 0): 


lim 

P.-KD 

1 ^1 

1/2 

l_i cSR 

R + L - 2S 

lim 

H“>00 

1 d|i _ 

1/2 


R + L - 2S 

and 

lim 

(i— *CO 

1 

1 

[i c5? ~ 

2SP 

2P _ _ 


So, at resonance (and ijj f 0) : 


_i 

M. 3n7 
j 
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If, on the other hand = 0: 


, 2 = l 


So 


1 ^ _ 1 I^£ = I^L = 0 

u i; 2L ’ m3r ll dP - 


and therefore: 


l 

M- . 


= - lim 


= - 11m 


L 2N .L 
J J 


L.-« 2N.|1 - L. - S L . \ 
J J Wj xj 


Consequently, if = 0 at resonance: 

1^1 1 _ 

(J. cJnT ~~ 2N . 


From the above limits and Eq. (C.8) it can be seen that, for any wave 
normal angle: 

2 > 2 


1 
M- 


-TTT =* ± CO (if R . < L . ) 
OB J J 


at resonance. 
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APPENDIX D. CURVES OF THE COMPONENTS OF dfr/ds FOR 
THE SAME CONDITIONS AS IN CHAPTER 2, 

BUT FOR i|f = -89° 

/ 

The following curves in Figs. 42 to 45 represent the several terms 

of d\ji/ds for both propagating modes in the same conditions as illustrated 

in Figs. 7 to 10 in Chapter 2, but for a large wave normal angle (\)t = -89°) 

The main differences from the small wave normal example examined in 

Chapter 2 are in the slow mode, due to the decrease of the resonance level 

o 

to an altitude below the crossover altitude. If i|t = 90 this resonance 
level would have been reduced to the two-ion resonance level also marked 
in the figures. In this example, the tendency of the particle density 
term toward infinity at resonance (see Fig. 42) is more apparent than 
in the case of small wave normal angle. \ 
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RAD/km 


FIGURE 42. CURVES OF THE COMPONENT OF di|(/ds DUE TO THE 
ELECTRON AND ION DENSITIES (N . ) , FOR A WAVE 
PROPAGATING IN THE SLOW MODE. 3 


SEL-71-055 


- 106 




0~ 4 10 -3 10 2 10 


RAD/km 


FIGURE 43. CURVES OF THE COMPONENTS OF d^/ds DUE TO 

THE MAGNITUDE (B) AND INCLINATION (p) OF THE 
EARTH'S MAGNETIC FIELD, FOR A WAVE PROPAGATING 
IN THE SLOW MODE. 
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RAD/km 


FIGURE 45. CURVES OF THE COMPONENTS OF dijr/ds DUE TO THE 

MAGNITUDE (B) AND INCLINATION (g) OF THE EARTH'S 
MAGNETIC FIELD, FOR A WAVE PROPAGATING IN THE 
FAST MODE. 
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